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Nuclear  Magnetic:  Resonance  in  Imperfect  Crystal*  * 

by 

N.  Bloembergen 

Division  of  Applied  Science,  Harvard  University 
Cambridge,  Massachusetts 


Abstract 


A systems** c survey  is  made  how  various  kinds  of  imperfections  in 
a crystalline  lattice  will  affect  the  position,  breadth,  shape  and  relaxation 
time  of  the  nuclear  magnetic  resonance.  In  accordance  with  Seitz! a classi- 
fication the  following  imperfections  are  considered:  (a)  dislocations; 

(b)  vacant  lattice  sites  and  interstitial  atoms  ; (c)  foreign  atoms  in  either 
interstitial  or  substitutional  position;  (d)  electrons  and  holes;  (e)  phonons; 
(f)  excitons.  Their  interaction  with^the  magnetic  dipole  moment  and  the 
electric  quadrupole  moment  of  the  nuclei  at  the  normal  lattice  sites  is  dis- 
cussed and  the  available  experimental  information  is  reviewed. 


1 

Introduction 


The  nature  of  imperfections  in  nearly  perfect  crystals  has  been  studied 
extensively  during  the  last  two  decades.  An  excellent  synthesis  of  this  broad 
field  has  been  given  by  Seitz  [ 1]  . At  the  time  his  report  was  written  only  one 
or  two  connections  between  crystalline  imperfections  and  nuclear  magnetic 
resonance  were  known.  During  the  past  few  years  several  new  examples  of 
the  interplay  of  magnetic  resonance  phenomena  «od  imperfections  in  a crys- 
talline lattice  have  been  found.  Some  of  these  have  recently  been  reviewed 
by  Pound  [2] , some  other's  will  be  discussed  in  more  detail  during  the  Bristol 
Conference . It  may  be  profitable  to  give  a broad  and  systematic  review  of 
the  possible  interactions  between  nuclear  spins  and  lattice  imperfections.  In 

♦ Part  of  this  research  was  supported  jointly  by  the  Office  of  Naval  Research, 
Army  Signal  Corps  and  the  U.3.  Air  Force.  This  paper  contains  the  material 
cf  two  lectures  presented  at  the  Bristol  Conference,  July  1954.  The  author  is 
indebted  to  the  ONR  for  transatlantic  transportation. 
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the  course  of  this  survey  the  opportunity  will  arise  from  time  to  time  to  point 
unsolved  problems  and  directions  in  which  further  investigation  might  be 
fruitful.  Some  special  attention  will  be  paid  to  metals  and  alloys,  where 
conduction  electrons  and  foreign  solute  atoms  cause  prominent  effects.  It 
is  hoped  that  this  review  will  stimulate  further  research  in  a field,  the  ex- 
ploration of  which  seems  to  be  still  in  its  infancy. 

According  to  Seitz  the  imperfections  to  be  considered  can  be  classified 

as: 

(a)  dislocations 

(b)  vacant  lattice  sites  and  interstitial  atoms 

(c)  foreign  atoms  in  either  aubstitional  or  interstitial  positions 

(d)  free  electrons  and  holes 

(e)  lattice  phonons 

(f)  excitons 

These  imperfections  may  change  the  distribution  of  the  internal  magnetic 
fields  and  the  gradients  of  the  crystalline  electric  field.  In  Section  II  the 
interaction  of  these  imperfections  and  the  nuclear  magnetic  dipole  moments 
is  discussed,  and  Section  ill  deals  with  the  interaction  with  the  nuclear  electric 
quadripole  moments. 

*■  The  direct  effect  of  transient  imperfections  in  the  form  of  light  quanta, 
charged  and  uncharged  material  radiations  is  negligible,  unless  one  wishes 
to  consider  the  radiofrequency  field  inducing  the  nuclear  spin  transitions  it- 
self as  photons  disturbing  the  lattice. 


II 

Magnetic  Dipole  Interactions 

In  diamagnetic  crystals  the  internal  magnetic  field  acting  on  a nuclear 
spin  consists  almost  entirely  of  the  contribution  from  neighboring  nuclear 
magnetic  moments.  The  random  spin  orientation  of  these  neighbors  gives  rise 
to  a line  broadening  of  the  magnetic  resonance  line,  whose  second  moment 
is  given  by  the  celebrated  formula  of  V an  Vlcck  13] 
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where  the  second  moment  is  defined  by 

a.' 


Sv)^”  = f (v-v  )2  g(v ) dv  , f g(v)dv  = 1, 


g(v)  is  the  normalized  shape^actor  of  the  absorption  line.  The  first  summation 
on  the  right-hand  side  of  equation  (1)  extends  over  all  identical  nuclear  spins  j, 
the  second  sum  ever  all  different  spins  f.  The  nuclear  gyromagnetic  constant 
V relates  the  resonance  frequency  with  the  magnetic  field 


2»v0=TH0  . (2) 

rQj  and  r^  are  iniernuciear  distances  and  o is  the  angle  between  the  radius 

vector  and  the  magnetic  field  H0.  The  magnetic  field  acting  on  the  nucleus 

will  he  comoosed  of  the  external  field  H and  the  internal  or  local  fields. 

o 

The  dipolar  interaction  produces  random  variations  from  the  mean  value  HQ, 
which  find  their  quantitative  expression  in  equation  (1). 

A.  Dislocations 

In  the  neighborhood  of  a dislocation  the  iniernuciear  distances  will  have 
changed  somewhat,  but  usually  less  than  one  percent  The  corresponding 
change  in  the  width  of  the  resonance  given  by  equation  (1)  is  therefore  negli- 
gible. However,  if  the  nuclei  under  consideration  have  a quadrupole  moment, 
d Ciidiigc  iu  uxe  Vdiuc  uf  SCCGnd  rsCdSut  will  occur , as  will  b*  crpisisc u 
later  when  the  quadrupole  effects  are  discussed. 

The  contribution  of  the  electron  orbitals  to  the  magnetization,  which 
is  of  course  zero  for  completely  filled  bands,  is  usually  quenched  by  the  crys- 
talline electric  field  even  for  incompletely  filled  shells.  The  small  change  in 
the  crystalline  fields  due  to  ths  strains  around  a dislocation  will  not  change 
this  situation.  In  particular,  a diamagnetic  substance  will  remain  diamagnetic 
under  deformation  and  the  contribution  of  the  electrons  to  the  internal  magnetic 
fiek’  at  the  position  of  the  nuclei  in  the  solid  remains  negligible. 

It  is  conceivable  that  an  indirect  effect  m3y  occur  in  the  presence  of 


T- 


» 

& 


TR199  -4- 

paramagnetic  ions.  The  Hamiltonian  containing  the  crystalline  field  split- 
ting and  the  spin-orbit  coupling  of  these  ions  may  change,  and  therefore  the 
relaxation  time  of  the  electron  spins  could  be  changed  in  a strained  lattice 
if  the  symmetry  of  the  crystalline  field  is  lowered  by  the  deformation.  This 
might  in  turn  change  the  relaxation  time  of  the  nuclear  spins  by  a mechanism 
which  will  be  discussed  in  paragraph  (C)-  As  the  influence  of  deformation  on 
the  electron  spins  in  paramagnetic  materials  has  received  little  attention  and 
* it  is  only  a side  issue  for  the  nuclear  resonance,  we  shall  hot  pursue  this 
matter  any  further. 

We  arrive  at  the  conclusion  that  dislocations  in  general  have  a negligi- 
ble influence  on  the  nuclear  resonance  in  the  absence  of  quadrupole  effects 

(l  = 

B . Vacancies  and  interstitial  atoms 

The  presence  of  vacancies  and  interstitial  atoms  will  cause  an  obvious 
change  in  the  summations  nf  equation  (1).  Since  their  concentration  will 
usually  not  exceed  1 atomic  per  cent,  the  second  moment  in  the  static  lattice 
remains  practically  unchanged. 

Nevertheless,  vacancies  and  interstitials  can  have  a pronounced  effect  on 
the  line  width  and  the  spin  lattice  relaxation  time.  At  sufficiently  high  temper- 
atures they  make  diffusion  possible  and  as  the  individual  nuclei  carry  out  dif- 
fusion jumps,  the  internal  magnetic  field  acting  on  the  nuclear  spins  will 
have  a time  dependence  with  a characteristic  time  given  by 
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is  a constant  of  the  order  of  the  atomic  vibration  frequency,  E c 

, is  the 
act 


Z =vdiff  = vat  *** 

where  v 

is  the  energy  for  formation  of  a vacancy  or  interstitial  atom  and  E 

activation  energy  necessary  to  move  such  a vacancy  or  interstitial  over  an 

atomic  distance.  Time -dependent  internal  fields  have  been  discussed  by 

Bloembergen,  Purcell  and  Pound  [*4]  foi  nuclear  relaxation  effects  in 

liquids.  The  same  treatment  applies  to  motion  in  solids. 

The  line  width  due  to  dipolar  interaction  between  randomly  moving 

dipoles  with  a correlation  time  T is  given  by  the  approximate  formula 

1/2 
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and  the  relaxation  time  by 
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In  general  there  will  not  be  a single  characteristic  time,  and  a more 
accurate  integration  over  a distribution  of  7*'s  is  appropriate.  A refined 
analysis  for  the  diffusional  motion  and  its  effect  on  the  nuclear  relaxation 
has  been  given  by  Torrey  [5] . The  effects  of  diffusional  motion  in  solids 
on  the  nuclear  spin  resonance  was  first  discovered  at  the  University  of 
Illinois  in  metallic  sodium  f6]  , hut  similar  effects  will  occur  in  any  solid 
if  the  diffusion  is  rapid  enough.  The  main  features  of  equations  (4)  and  (5)  are 

(1)  When  the  diffusion  jump  frequency  - becomes  larger  than 

the  line  width  ^(6v)-^  ^ ^ of  the  rigid  lattice,  a motional  narrowing 
of  the  line  sets  in. 

(2)  The  relaxation  time  T,  will  pass  through  a minimum  when  v 
is  of  the  order  of  the  Larmor  precession  frequency  vq  of  the 
nuclear  spins.  By  measuring  the  line  width  and  T,  as  a function 
of  temperature  the  activation  energy  for  self-diffusion  is  deter- 
mined. The  diffusion  of  different  isotopes  can  be  studied  separately. 

A particularly  interesting  application  could  be  made  to  ionic  crystals, 
e.g.,  the  alkali  and  sil/er 'baUdes.  The  time  dependence  of  the  dipolar  inter- 
action between  two  slow-moving  anion-nuclei  would  he  determined  by  the  jump 
frequency  of  the  anions,  but  the  dipolar  interaction  between  cation-anion  and 
cation-cation  pairs  would  have  a frequency  spectrum  determined  largely  by  the 
Hgher  jarrp  frequency  of  the  cations.  For  the  nuclear  resonance  of  the  anions 
one  would  have  two  terms  on  the  right  side  of  equations  {4)  and  (5),  one  with 
a longer  T for  the  dipolar  interaction  between  two  anions,  and  one  with  a 
shorter  T for  the  anion-cation  dipolar  interaction.  Diffusion  effects  in  solids 
are  discussed  in  detail  in  papers  by  Andrew  and  Slichter  at  the  Bristol  Con- 
ference. 

An  anion-vacancy  occupied  by  an  electron  (F-cent er)  and  other  com- 
binations of  vacancies  and  other  imperfections  will  be  discussed  in  the  next 
paragraph. 

C.  Foreign  atoms 

A nonmagnetic  foreign  atom,  either  in  a regular  or  interstitial  site,  will 
have  a negligible  effect  on  the  magnetic  field  produced  at  the  position  of  sur- 
rounding nuclei. 
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Magnetic  impurities  will  produce  a large  local  field,  of  the  order  of 
3 

10  oersted,  at  nuclei  which  are  nearest  neighbors,  but  this  field  will  drop 
off  as  the  inverse  cube  of  the  distance.  Consequently  the  resonance  absorp- 
tion of  nuclei  in  the  immediate  neighborhood  will  be  displaced  and  contribute 
to  the  far  wings  of  the  absorption  line,  which  in  general  will  remain  unobserved. 
The  bulk  of  the  nuclei  will  be  practically  unaffected  and  the  second  moment  of 
the  observed  signal  will  remain  almost  the  same. 

Furthermore,  the  magnetic  field  of  the  impurity  is  not  static  but  changes 
with  the  spin  orientation  of  the  magnetic  ion.  The  characteristic  time  is  the 
electron  spin  lattice  relaxation  time  of  the  impurity.  For  paramagnetic  salts 
with  high  concentrations  of  magnetic  ions  the  characteristic  time  could  alter- 
nately be  determined  by  the  spin- spin  or  exchange  interaction,  whichever  is 
the  fastest.  For  isolated  magnetic  imperfections  the  spin-orbit  lattice  coupling 
is  the  only  mechanism  causing  variation  of  spin  orientation  with  time.  It 
has  been  shown  conclusively  [2,7,8]  that  in  many  ionic  and  valence  crystals 
the  nuclear  spin  lattice  relaxation  time  is  determined  by  the  presence  of 
these  magnetic  imperfections.  They  can  be  chemical  impurities  of  the  tran- 
sition group  elements,  F-centers,  filled  donor  or  acceptor  levels  in  semi- 
conductors, or  other  imperfections  provided  they  are  carrier*  of  electronic 
magnetism. 

The  nuclear  spin  relaxation  takes  place  by  the  following  sequence  of 
events  [7]  . An  absorbed  quantum,  h vQ*is  transported  by  spin  diffusion  through 
the  nuclear  spin  system  toward  the  neighborhood  of  a paramagnetic  impurity . 
This  process  is  independent  of  the  lattice  temperature.  The  quantum  hvo 
is  then  transferred  to  the  lattice.  The  lattice  vibrations  modulate  the  orbit 
of  the  magnetic  imperfections  and  via  the  spin  orbit  coupling  flip  the  elec- 
tronic spin  The  time -dependent  local  field  of  this  electron  spin  in  turn  flips 
the  nuclear  spin.  This  process  is  determined  by  the  electron  spin  lattice 
relaxation  time,  which  depends  on  the  lattice  temperature.  In  principle  the 
measurement  of  the  nuclear  relaxation  tin:  i as  a function  of  temperature  will 
give  information  about  the  concentration  of  paramagnetic  imperfre lions  and 
their  spin  lattice  coupling.  It  should  be  pointed  out  that  much  more  detailed 
information  is  obtained  from  the  observation  of  the  electron  spin  resonance 
of  the  impurity  itself.  The  significance  in  the  case  of  nucle  \r  resonance  is  . 
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Fig.  2.  A lattice  of  nuclear  spins  with  a paramagnetic  impurity.  The  reso- 
nance is  spread  out  over  the  range  of  local  fields  produced  by  neighboring  di- 
poles. A quantum  hvQ,  absorbed  by  spin  A,  is  transferred  to  the  site  B by  a 
simultaneous  spin  flip.  The  magnetic  energy  quantum  diffuses  toward  the  im- 
purity. It  is  transferred  to  the  lattice  via  interaction  with  the  electron  spinS 
and  the  spin -orbit  coupling. 
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mainly  that  these  impurities  provide  an  adequate  relaxation  mechanism. 


D.  Electrons  and  holes 


Since  free  electrons  and  holes  are  carriers  of  an  electronic  spin,  they 
create  internal  magnetic  fields  which  interact  with  the  nuclear  spins.  It  is 
a matter  of  taste  whether  one  wishes  to  consider  the  conduction  electrons  in 
a metal  as  imperfections , but  in  semiconductors  this  is  certainly  the  generally 
accepted  point  of  view.  To  achieve  a uniform  presentation  conduction  electrons 
will  be  considered  as  imperfections  in  this  paper. 

In  a metal  with  Fermi  energy  Eg.  most  of  the  electron  spins  are  paired 
off  and  only  a fraction  of  the  order  of  SH^/rX.  has  unbalanced  magnetic  mo- 
meits.  The  Pauli  paramagnetism  oi  the  conduction  electrons  gives  rise  to 
a shift  of  the  nuclear  resonance  in  metals  to  higher  frequencies,  first  dis- 
covered by  Knight  [9]  . The  relative  shift  is  given  by  [10,11] 
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»b?r*  5 is  the  Bohr  t> — A I - S is  the  hypexfine  interaction  in  the  cor- 
responding atomic  state; ^ is  the  square  of  the  electronic  wave 
function  in  the  metal  evaluated  at  the  position  of  the  nucleus,  normalised  over 
the  atomic  volume  and  averaged  over  the  Fermi  surface.  N(Ep.}  is  the 
density  of  states  of  the  conduction  electrons  at  the  Fermi  surface,  per  unit 

volume  and  per  unit  energy  interval.  The  observed  relative  shifts  range 

7 joo 

from  0.0261  percent  for  Li  to  2.5  percent  for  Hg  ' . There  is  a general 
increase  with  increasing  atomic  number,  as  the  hyperfine  interaction  becomes 
larger. 

In  alloys  the  shift  given  by  equation  (6)  ma;  vary  [11]  and  will  indicate 
a variation  in  the  product  of  the  density  of  states,  the  atomic  volume  and  the 
probability  density  of  the  electrons  at  the  nucleus.  In  fact,  the  pr  cdi_t  may 
show  a dispersion  and  have  different  values  at  different  nuclei  depending  on 
the  relative  position  of  the  impurities  in  the  neighborhood.  In  general  this 
will  give  rise  to  a spurious  asymmetric  broadening  of  the  line. 

Mctt[  12}  and  Friedel  [13]  have  calculated  the  shielding  of  the  extra 
charge  on  the  Zn  atom  oy  the  cuuCcntisuun  of  conduction  electrons  around 
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it.  The  ihieWing  i»  practically  complete  already  at  the  lmi**?  neighbors 
and  since  the  impurity  was  introduced  as  a neutral  atom,  tbs  electron  density 
throughout  the  bulk  of  the  metal  ia  unchanged. 


However,  the  Knight  shift  is  not  determined  by  the  total  density  of 
conduction  electrons  bat  by  j 4>  ^{5}  evaluated  at  the  Fermi  surface.  One- 
electron  wave  functions  with  a fixed  magnitude  of  the  wave  number  k must 
be  evaluated  in  the  neighborhood  of  the  impurity . Brooks  * has  pointed  out 
that  in  the  same  free  electron  approximation  as  used  by  Mott  the  square  of 
the  wave  function  for  electrons  with  a given  value  of  k falls  off  more  slowly 
than  the  exponential  Law.  The  value  k corresponding  to  the  — sn-^hnri  a* 
the  Fermi  level  must  be  taken.  Matt's  result  of  exponential  shielding  is  ob- 
tained by  integration  over  all  occupied  values  of  k.  Reliable  values  for  the 


shift  can  only  be  obtained  if  the  wave  equation  in  the  periodic  crystalline 


yuicutial,  ycltiuvcu  W|  the  liujiiujiiy,  la  Svtvcd.  It  IS  dca-i  that 
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will  be  a function  of  r,  the  distance  from  the  impurity.  It  may  show  some 
oscillations  before  leveling  off  to  its  value  in  the  bulk  of  the  metal. 
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For  very  low  impurity  concentrat  ions  the  main  absorption  peak  should 
remain  undisplaced  with  tome  small  {usually  unobservable.)  _ attributions  on 
either  side  from  nuclei  near  the  i mp" t > iy , For  concentrations  larger  than 
a few  percent,  a broadening  of  the  line  is  caused  by  the  distribution  of  shifts, 
and  the  frequency  of  maximum  absorption  may  change  since  the  majority  of 
the  nuclei  are  now  near  an  impurity.  These  effects  have  been  observed 
qualitatively  for  solutions  of  tin  in  thallium,  as  shown  in  Fig.  4.  There  is, 
however,  some  uncertainty  about  the  phase  diagram  and  a two-phase  region 
may  exist  between  4 and  10  percent  atomic  percent  tin.  Furthermore,  the 
line  broadening  wa*  not  asymmetric  as  might  in  general  be  expected.  If 
the  increase  in  width  is  indeed  due  to  a dispersion  of  Knight  shifts,  and  there- 
fore proportional  to  vft,  the  increase  will  be  more  pronounced  in  high  fields. 
More  precise  measurements  on  various  alloy  systems  with  spin  1 — 1/2  are 
badly  nee  dew. 


it  is  a matter  of  taste  whether  one  wishes  to  consider  the  alloys  with  a 
high  concentration  a?  solute  atoms  as  imperfect  crystals.  As  there  is  no  clear- 
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cut  concentration  at  which  an  impurity  becomes  a regular  component  of  the 
solid,  alloys  will  be  considered  as  imperfect  solids  in  this  review  to  allow 
for  a convenient  general  presentation  of  the  subject. 

Whereas  the  Knight  shift  changes  continuously  with  composition  in  a 
solid  solution*  discontinuous  changes  are  expected  in  going  from  one  phase 
to  another.  In  two-phase  regions  the  nuclear  resonance  is  doubled.  This 
has  been  observed  for  the  Tt205  resonance  in  the  thallium-mercury  system 

[H]. 

In  the  liquid  phase,  and  perhaps  also  just  below  the  melting  point, 
the  nuclei  will  diffuse  rapidly  and  the  resonance  frequency  will  be  deter- 
mined by  the  average  value  of  the  product  vj  »{*__,.(  0)  j 2 N(E_,),  where  the 
average  is  taken  over  all  configurations  around  the  nucleus  in  question. 

Instead  of  the  dispersion  of  Knight  shifts,  the  center  of  gravity  is  measured 
in  the  liquid.  An  asymmetric  impurity  broadening  in  the  solid  corresponds 
to  a variation  of  the  resonance  frequency  in  liquid.  This  last  effect  is,  in 
terms  of  Seitz' s classification,  the  combined  influence  of  three  types  of 
imperfections:  electrons,  foreign  atoms  and  vacancies. 

The  shift  caused  by  conduction  electrons  in  semiconductors  is  negligible 
[14].  The  total  density  of  electrons  and  also  of  those  with  unbalanced  spin 
is  so  much  smaller  than  in  metals  that  no  observable  effect  results. 

The  nuclear  relaxation  in  metals  iB  usually  also  determined  by  the 
interaction  of  nuclear  spins  and  conduction  electron  spins.  A simple  picture 
is  that  the  nuclear  spin  has  a certain  probability  to  make  a transition,  when- 
ever =n  electron  in  the  electron  gas  passes  by.  Only  a small  fraction  cf  the 
electrons,  of  the  order  of  kT/E^,,  representing  the  tail  of  the  Fermi  distri- 
bution, are  allowed  to  take  x.p  the  small  energy  differences  involved  in  the 
magnetic  transition.  The  relaxation  time  should  be  independent  of  the 
external  field  HQ,  although  experimentally  some  variations  of  T^  with  Hq 
have  been  observed.  If  one  lakes  into  account  only  the  point -interaction 
caused  by  the  finite  value  of  the  wave  function  at  the  position  of  the  nucleus 
the  following  relation  due  to  Korringa  [15]  exists  between  the  relaxation  ime 
and  the  shift. 

/Av.  \Z  „ 2 
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The  relaxation  time  gives  the  product  of  the  same  quantitie  s va| 
as  the  shift. 

Actually  the  part  of  the  electronic  wave  function  with  non  S-character 
{p-,  d-,  etc. ) will  also  contribute  to  the  relaxation  time,  but  not  to  the  shift. 

The  equal  sign  in  equation  (7)  should  therefore  be  replaced  by  a<C  sign.  This 
correction  of  the  Korringa  relation  is  probably  more  important  than  the  neglect 
of  correlation  and  exchange  effects  in  the  electron-gas  model.  The  observed 
deviations  from  the  Korringa  relation,  which  are  mostly  in  the  right  direction 
[16]  , could  perhaps  find  an  explanation  in  this  fashion. 

Although  the  hyperfine  interaction  of  p-  and  d- orbits  is  usually  much 
smaller  than  in  S-states,  the  amount  of  S-character  of  the  one-electron  wave 
functions  may  be  small  even  in  alkali  metals.  Jones  and  5c niff  [17]  estimate 
that  in  sodium  there  is  60  per  cent  S-character  and  in  lithium  25  per  cent. 
Assuming  the  hyperfine  interaction  in  the  p-state  to  be  only  10  per  cent  of  that 
in  the  S-state,  the  interaction  of  the  p-part  would  reduce  the  relaxation  time 
calculated  from  the  Korringa  relation  by  4Qper  cent  in  lithium.  A careful 
systematic  check  of  this  relation  is  desirable. 

The  relaxation  time  in  semiconductors  will  usually  be  determined 
by  electron  spins.  It  is  predicted  that  the  spins  of  bound  electrons  or  holes 
which  act  like  paramagnetic  impurity  atoms  will  produce  a more  efficient  re- 
laxation mechanism,  discussed  in  the  preceding  paragraph,  than  the  spins  of 
free  carriers  [ 14]  . No  experimental  evidence  is  available,  but  in  principle 
the  nuclear  relaxation  could  give  a completely  independent  determination  of 
the  number  of  bound  electrons  or  free  carriers  in  semiconductors. 

In  noncubic  crystals  the  p-,  d-  etc.  character  of  the  electron  orbitals 
gives  risk  to  an  anisotropy  of  the  Knight  shift,  which  has  been  observed  [11] 
in  white  tin  and  solid  mercury.  Figure  6 illustrates  how  an  ellipsoidal  dis- 
tribution in  space  of  the  density  of  electrons  with  unbalanced  spins  will  give 
rise  to  an  aiding  or  counteracting  local  field  on  the  nucleus,  depending  on  the 
direction  of  Hq  with  respect  to  the  crystallographic  axes  (=  axes  of  the  electron 
density).  In  polycrystalline  samples,  which  must  necessarily  be  used  in 
resonance  work  on  conductors,  the  line  gets  a characteristic  asymmetric 
shape.  The  observed  anisotropy  is  proportional  to  the  field  strength,  as  it 
increases  proportionally  to  the  number  of  electrons  with  unbalanced  spin  in. 
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the same  manner  as  the  isotropic  shift. 

*V  1 * n<  V ^iB  (8) 

o sav 

The  first  term  arises  from  the  deviation  from  cubic  symmetry 

qF=  f*Fr~5(3z2  - t2}**  dV 

of  the;  density  of  electrons  near  the  Fermi  level.  The  second  term  is  caused 

by  the  anisotropy  in  the  g-factor  of  the  conduction  electrons.  Unfortunately 

we  have  only  one  relation  between  q^,  and  g||  - g^.  Both  quantities  are 

difficult  to  determine  theoretically  and  experimentally  by  other  means.  If 

the  anisotropy  is  small,  the  effects  can  be  made  observable  in  very  high 

fields  H . 

o 

In  alloys  one  will  have,  in  addition  to  the  distribution  of  isotropic  shifts, 
a distribution  of  anisotropic  shifts.  This  will  add  to  the  broadening  of  the 
resonance  In  such  systems. 

' Although  there  is  no  experimental  evidence,  we  wish  W point  out  the 
possibility  of  investigating  order -disorder  transitions  in  alloys  in  this  manner. 
In  the  first  place  there  will  be  an  obvious  change  in  the  dipolar  contribution  to 
the  line  width  on  ordering,  which  can  be  calculated  with  equation  (1).  This 
change  in  width  will  be  especially  marked  if  the  other  constituent  in  the  binary 
alloy  has  zero  spin.  In  the  second  place  the  dispersion  in  the  Knight  shift  will 
cease  to  exist  in  the  completely  ordered  state.  Ordering  will  in  general  be 
accompanied  by  a marked  change  in  the  lirf.  width. 

Since  the  width  of  the  nuclear  resonance  is  largely  determined  by  the 
near  neighbors,  it  is  a good  indicator  for  short-range  order  without  the 
complications  of  temperature  diffuse  scattering  which  are  inherent  in  the  X-ray 
technique. 

The  onset  of  precipitation  will  also  cause  a change  in  composition  of 
the  immediate  environment.  The  precipitation  hardening  in  a system  with 
solvent  atoms  with  1 = 0 and  solute  atoms  with  1=  1/2  may  be  studied  profitably. 
The  nuclear  resonance  will  broaden,  when  a phase  rich  in  solute  atoms  precipiti- 
tate  s- 
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Lattice  vibrations  will  modulate  the  dipolar  interaction  between  nuclear 
spins  and  thus  produce  a relaxation  mechanism.  In  fact,  this  was  historically 
the  first  magnetic  spin  lattice  relaxation  mechanism,  proposed  by  Waller  [17] 
in  1932  for  electron  spin  systems.  It  can  be  shown  that  nuclear  relaxation 
times  calculated  by  this  theory  are  in  general  a factor  of  10^  to  10^  longer 
than  observed.  The  dominant  relaxation  mechanism  in  metals  is  caused  by 
conduction  electrons;  in  ionic  and  valence  crystals  frequently  by  paramagnetic 
imperfections.  The'  energy  in  the  former  case  is  taken  up  by  a change  in 
momentum  of  the  electrons;  in  the  latter  case  the  lattice  vibrations  are 
essential  to  the  transfer  of  energy.  The  relaxation  is  thus  strictly  -peaking 

- a _ rr ,j.  ^ c : £ i.: 

CL  LUlUUIilCU  CUCL^  U1  tWW  o. 

In  many  molecular  and  certain  other  crystals  different  types  of  thermal 
motion  than  lattice  vibrations  occur,  such  aB  hindered  rotation  tunneling 
between  two  or  more  positions  of  equilibrium.  Such  motions, which  produce  a 
very  effective  relaxation  mechanism  for  nuclear  spins,  should  perhaps  also 
be  considered  as  imperfections  of  the  perfect  rigid  lattice.  A wealth  of  ex- 
perimental data  exists  on  this  matter,  for  which  the  reader  iB  referred  to 
the  literature  [18]  . In  this  context  the  diffusion  phenomena  discussed  under 
(b)  should  be  considered  as  the  combined  effect  of\Eacancies  and  phonons. 

F.  Excitons 

An  excited  atom  or  ion  in  a solid  will  have  no  effect,  if  the  excited 
state  is  nonmagnetic.  If  the  state  is  magnetic,  it  will  have  a similar  effect 
on  nuclear  relaxation  as  a paramagnetic  impurity.  The  characteristic  time 
of  its  local  field  will  be  determined  by  the  life  time  of  the  excited  state,  by 
its  diffusion  through  the  lattice  or  by  itB  spin  reversal  time,  whichever  process 
is  the  fastest.  No  experimental  evidence  is  known  to  the  author.  The  con- 
centration of  excitons  is  usually  too  low,  but  a change  In  nuclear  relaxation 
time  in  phosphors  under  strong  illumination  might  be  detectable. 

Ill 

Quadrupole  Effects 

If  the  nuclear  spin  I — 1,  the  interaction  of  the  imperfection  with  the 
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Fig.  7.  The  observed  asymmetric  nuclear  resonance  in  polvcry sta Tine 
p~tin.  The  anisotropy  in  the  Knight  shift  is  proportional  to  the  resonant 
frequency.  The  asymmetry  is  more  pronounced  in  the  recorded  derivatives . 
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electric  quad ru~ ole  moment  will  frequently  be  more  important.  The 
crystalline  electric  field  will  always  be  disturbed  by  the  imperfection,  and 
special  magnetic  properties  of  unbalanced  electron  spins  are  not  required. 

The  elements  of  quadrupole  interaction  m an  inhomogeneous  electric 
field  of  axial  symmetry  will  first  be  recalled  [19>  20] . Assume  that  the  quad- 
rupole interaction  can  be  considered  as  a small  perturbation  on  the  21  + 1 
equally  spaced  energy  levels  in  the  magnetic  field  Ho>  which  is  taken  in  the 
z-direction.  The  equal  xpacing  will  be  destroyed.  Introduce  a scalar 
quantity 

eq  = y~e.(3cos2d  ■ -l)r‘3 

j 


where  dj  is  the  angle  between  the  radius  vector  rj  toward  a charg  c 6 j 1X1  tuC 
lattice  and  the  axis  of  symmetry  at  the  origin.  Denote  the  angle  between  this 
symmetry  axis  and  the  magnetic  field  by  8 . Then  the  energy  difference  be- 
tween the  levels  nij  and  -1  of  the  nuclear  spin  is  given  by  first-order 
perturbation  theory. 


A £ , = h v +(2m  - 1)(3  cos  6 - 1) Avf**  ^ 5?  (?) 

m — - 1 O • i; 

where  Q is  the  nuclear  quadrupole  moment.  It  is  seen  from  c ^nation  (9) 
that  the  distance  between  the  levels  m^  = 1/2  and  - 1/2  is  unperturbed  in  this 
approximation.  Second-order  perturbation  theory  gives  for  the  energy 
difference  between  these  levels  if  the  nuclear  spin  ie  an  odd  half-integer. 


^1/2-*-  1/2 


hvo 


21  + 3. 
4IZ(2I-1) 


402  2 

e „Q 

hv 


(1-9  cos2  0 )(1  - cos20  ) 


(10) 


The  splitting  of  the  nuclear  resonance  has  been  observed  by  Pound  [19]  and 
others  [20]  in  noncubic  single  crystals.  In  cubic  crystals,  or,  more  accu- 
rately, at  points  around  which  the  symmetry  i»  cubic,  q - G and  there  is  ro 
quadrupolar  interaction.  In  poly  crystalline  noncubic  crystals  the  resonance 

of  the  satellite  lines,  i.  e.»  the  transitions  excluding  the  central  component 

2 

mj  = 1/2— >-  1/2,  will  be  spread  over  a large  frequency  interval,  as  3 cob  0-1 
ranges  from  it«  maximum  value  of  2 to  its  minimum  value  of  - 1.  An  inte- 
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gration  over  the  solid  angle  sinO  dO  d^'  yields  for  the  absorption  in  a 
powdered  sample  the  line  shape  [21]  shown  in  Fig.  9.  In  second  approxi- 
mation even  the  central  component  will  be  broadened  and  its  asymmetric 
shape  for  a poly  crystalline  sample  is  given  in  Fig.  10.  On  these  curves 
the  dipolar  broadening  ought  to  be  superimposed  . A similar  curve  had 
been  given  by  Nierenberg  and  Ramsey  [ 22]  in  molecular  beam  experiments 
on  alkalihalide  molecules.  The  observation  technique  of  nuclear  magnetic 
resonance  makes  it  very  difficult  to  detect  the  broad  low  tails  of  these  pat- 
terns. With  the  usual  modulation  method  the  slope  of  the  absorption 
curve  is  detected  and  this  quantity  varies  inversely  as  the  square  of 
broadening.*  Frequently  only  the  absorption  near  the  resonance  peak  can 
be  observed.  The  quantity  <X(vrnax)  will  decrease  as  the  quadrupole 
interaction  increases,  the  total  area  under  the  absorption  curve  remaining 
constant.  It  will  be  useful  to  introduce  normalized  shape  functions  lot  - 
each  satellite  and  for  the  central  component  separately.  For  the  complete 
line  we  have 

- ft  a » /..\  x ft  £ » f..v  r — .t-3  ;in 

o\ ' i -•  * erv»(  • * ■ ” s *"»  J y \ * */ 


*,v>  = 9/35  gc(v|  +r4  g3/2,l/2(v)  + X g5/2,3/2(v) 

3 35-8-1/2,-3/2(v)  + ~3T  g-3/2r5/2lv) 
for  I - 5/2,  etc. 

In  crystals  with  noncubic  symmetry  a splitting  of  the  nuclear  spin  levels 
occurs  even  in  the  absence  of  a magnetic  field,  as  indicated  in  Fig.  8. 
Transitions  between  these  levels. give  rise  to  so-called  pure  quadrupole 
resonance  line*  [ 23],  although  it  should  be  remembered  that  the  transi- 
tions are  always  induced  by  a magnetic  radiofrequency  field.  A single 
crystal  is  not  needed,  as  the  magnetic  field  is  absent.  In  powders  the 
lines  will  broaden  immediately  upon  application  of  an  external  field  We 
shall  now  discuss  the  influence  of  the  various  imperfections. 


♦ When  a square -wave  modulation,  on  -and  off-resonance,  is  used,  the  signal 
is  inversely  proportional  to  the  broadening  itself. 
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Fig.  8.  Energy  levels  of  a nuclear  spin  with  X = 3/2.  in  a magnetic  field  with 
cubic  environment,  first-order  perturbation  in.  an  axially  symmetric  field  and 
the  zero  magnetic  field  splitting. 


Fig.  9.  First-order  quadrupole  perturbation.  Line  shape  for  I = 3/2  in 
powdered  samples  of  axially  symmetric  crystals  (dcished  line).  With  dipolar 
broadening  superimposed  the  drawn  curve  results.  Frequently  th2  satellites 
are  spread  out  ove:r  such  a large  frequency  range  that  the  wings  become  un- 
observable. 
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A , Dislcctitiors 

The  strain  field  around  a dislocation  line  will  produce  a change  in 
the  gradient  of  the  electric  field  at  the  nuclei  near  ?he  dislocation.  The 
change  m gradient  Aq  will  be  proportional  to  the  strain,  which  in  turn  is 
inversely  proportional  to  the  distance  from  the  dislocation  line  A distribu- 
tion of  gradients  q and  of  angles  8 will  be  created.  Actually  the  nuclei  will 
be  at  points  where  the  axial  symmetry  is  lost.  Formulae  for  the  quadrupole 
interaction  in  fields  of  lower  symmetry  have  been  worked  out  [20]  . Qualitatively 
the  results  for  a distribution  of  gradients  with  or  without  axial  symmetry  will 
be  the  same.  A general  broadening  of  the  components  of  the  resonance  will 
occur.  The  central  component  will  only  be  affected  in  second-order  approxi- 
mation. In  noncubic  single  crystals  dislocations  will  broaden  the  satellites 
much  more  than  the  central  component.  The  latter,  which  is  of  course  only 
present  if  the  nuclear  spin  is  an  odd  multiple  of  one  half,  is  broadened  in- 
versely proportional  to  the  magnetic  field  Hq, 

The  strains  will  distort  the  euhie  symmetry  in  cubic  crystals.  Before 
the  deformation  all  components  coincide  in  the  single  crystal  or  polycry  stalling 
sample  of  cubic  symmetry.  After  the  deformation  the  satellite  components 
will  be  broadened  so  much  that  they  are  frequently  unobservable , ard  only  the 
central  component  remains.  Watkins  [24]  discovered  this  effect  for  the 
bromine  and  iodine  resonances  in  rather  perfect  crystals  of  KJBr  and  KI.  He 
introduced  a calibrator  into  the  experimental  apparatus  which  ma^e  (absolute) 
intensity  measurements  possible.  The  intensity  of  these  resonances  cor- 
responds to  that  of  the  central  component  alone,  namely  0.4  of  the  total  in- 
tensity  for  the  bromine  isotopes  with  spin  I = -j,  and  yg-  for  iodine  with  I = -j. 
When  the  crystals  were  subjected  to  severe  cold  work  (plastic  flow  resulting 
in  a 22  per  cent  charge  in  dimension  along  the  100  axis)  the  asymmetric 
broadening  of  the  central  component  was  alpo  observed.  Superimposed  on 
the  plastic  deformation  effect  was  a reversible  broadening  of  3 per  cent  to 
elastic  strains.  The  order  of  magnitude  of  the  quadrupole  interaction  in  the 
stress  field  of  these  dislocations  may  be  obtained  from  the  ullowing,  admit- 
tedly very  crude,  argument. 

Consider  first  the  change  in  electilc  tieid  gradient,  if  a neighboring 
charge  e at  distance  r is  displaced  over  a distance  Ar.  One  finds  for  the 
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the  change  in  gradient.due  to  this  compressional  strain  c = - Ar/r 

■A  (eq)  = 6e  r \ 

The  actual  gradient  at  the  nucleus  is  not  produced  by  the  displacements  of 
bare  charges.  The  gradient  at  the  nucleus  is  a very  sensitive  function  of 
the  distortion  of  the  ion  core  around  it.  Relatively  slight  distortions  of  this 
core  can  produce  large  distortions  in  the  symmetry  of  the  wave  function  in 
the  immediate  neighborhood  of  the  nucleus.  Since  the  electron  density  is 

_3 

weighted  by  a factor  r in  the  nuclear  quadrupole  interaction,  the  actual 
gradient  at  the  nucleus  may  be  increased  by  an  antishielding  factor  l-f|y 
which  has  been  calculated  by  Foley  and  Sternheimer  [?,5]  for  a number  of 
ions.  Tjhey  find  that  this  factor  1 x!y_!  is  47.5  for  d-  , 9.7  for  Gu+. 

51.2  for  Rb  and  98.8  for  Cs  . The  order  of  magnitude  of  the  quadrupole 
interaction  is  changed  by  this  factor 

The  ion  core  around  a given  nucleus  can  be  distorted  not  only  by  the 
displacement  of  external  charges,  but  also  by  the  interaction  with  neigh- 
boring ion  core 8.  In  general  there  will  be  an  interplay  between  the  charge 
or  "valance11  effect  and  the  "sine"  effect  ef  the  ion  cores.  Both  effects  will 
presumably  be  important  in  the  case  of  crystal  deformation.  A generalized 
multiplication  factor  X is  introduced  for  each  crystal  lattice,  which  takes 
account  of  the  ion  core  distortion  in  the  strained  lattice,  And  we  write  for 
the  change  in  quadrupole  interaction 

A (eq)  = 6Xe  a *Aa  = 6Xe  a ^ t 

where  a is  the  distance  between  nearest  neighbors  and  j is  a measure  cf  the 
strain.  The  effect  of  all  neighbors  is  incorporated  in  X,  which  is  a dimen- 
sionless constant  (more  precisely  a fourth-order  tensor)  in  the  assumed 
proportionality  between  the  strain  and  the  change  in  gradient  at  the  nucleus. 

For  a cubic  lattice  under  hydrostatic  pressure  the  net  effect  would  of  course 
be  zero,  but  at  a distance  r from  a dislocation  line  with  Burgers?  vector  b 

_ i 

we  may  estimate  the  order  of  magnitude  of  * as  br  The  stress  field 
around  tho  dislocation  lijoe  drbpn  off  inversely  proportional  to  the  distance  [26]  . 
As  b is  of  the  same  magnitude  as  a,  we  obtain 

v -2  -1 

* (eq)  = 6X  e a r 
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The  corresponding  change  in  resonant  frequency  of  the  transition  between  the 
nuclear  spin  levels  and  m^  - 1 is  consequently  in  first-order  perturbation 

hAv  = + (2m-  l)(3cos2  e-1)  -18?.Q-\—  (12) 

81(21  - \)tslC 

A good  measure  for  the  magnitude  of  the  perturbation  is  obtained  by 
calculating  its  second  moment  aB  a radial  average  over  the  deformed  cylinder 
around  a dislocation  line.  The  lower  limit  of  integration  is  rather  arbitrarily 
chosen  as  r.  = 3a.  The  second  moment  will  depend  only  logarithmically  on  this 
cut-off  radius.  Inside  this  radius  the  atomic  displacements  cannot  be  calculated 
from  a continuum  elastic  theory.  Moreover  the  relative  fraction  of  nuclei  in.- 
aide  this  radius  is  small,  and  they  c sntri  bute  to  the  far  wings  of  the  absorption 
line.  Their  signal  will  escape  observation  and  a cut-off  is  dictated  by  this 
physical  interpretation.  Th*  outer  radius  1-  chosen  so  that  the  cylinder  has 
the  average  volume  per  dislocation  line,  irr2  = c"1, where  c is  the  number  of 
dislocation  lines,  crossing  1 cm'  of  normal  area.  The  second  moment  of  the 
line  is  then  found  to  be 
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This  result  differs  by  a factor  3/ir  from  a more  rigorous  formula  derived 
by  Watkins  for  the  effect  of  six  neighboring  ions  due  to  shear  strain  in  ine 
K Br  type  lattice. 


(5V)2  = r_3efQ^2m-12|  x 2 c a2  ln  ( ! / 9 ca2)  ( 13) 

5ir  L 8I(2I-l)a  hJ 

random 

disl. 

Watkins  took  proper  account  of  the  angular  dependence  of  the  strain  a, round 
the  dislocation  line  [26]  and,  averaged  over  a statistical  distribution  of  dis- 
locations with  random  orientation,  ca  is  the  probability  that  a dislocation 
line  passes  through  a given  lattice  point.  Present  uncertainty  in  the  anti- 
shielding factor  X and  the  experimental  data  hardly  justify  a acre  detailed 
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calculation,  and  Watkim's  formula  may  be  used  equally  well  for  edge  and 
screw  dislocations.  The  phenomcn.oziolcgics.1  factor  X is  a function  of  the  ion 
and  the  lattice  in  which  it  is  embedded.  Its  order  of  magnitude  is  not  neces- 
sarily the  same  as  the  Stexnusuner  fa.'to r 1 t j ~y  { which  describes  the  dis  = 
tortion  of  an  ion  core  by  an  external  charge  and  is  a constant  for  a given  ion. 

Watkins  found  that  his  results  were  consistent  with  X = 10  for  the  bromine 

8 -2 

and  iodine  resonance  in  KBr  and  KI,  with  c = 6.  1C  cm  for  the  pure  crys- 
10  -2 

tals  and  c = 10  cm  for  the  plastically  deformed  crystals.  It  would  be  oi 
interest  to  compare  the  values  of  X for  anions  and  cations  in  the  same  crystal 
e.g.,  in  Na  Br.  It  is  to  be  expected  that  X should  be.  considerably  smaller 
for  Na+  than  for  Br  . 

An  effect  of  dislocations  on  the  satellites  has  been  found  in  metallic 
copper  [ 11].  The  specimen  consisted  of  320  mesh  copper  filings  to  permit 
the  radiofrequency  field  to  penetrate  the  individual  particles.  When  cold- 
worked  filings  were  annealed  the  peak  intensity  ot  the  resonance  could  be 
increased  by  a factor  2.5.  By  annea&xig  8s(vTOax)  in  equation  (11)  changes 
from  essentially  zero  to  gr(vrnax)»  which  regains  the  same.  For  intermediate 
stages  of  anneal  intermediate  peak  intensities  g(vmax)  were  recorded.  In 
the  well-annealed  crystals  the  density  of  dislocations  is  so  low  that  first-order 
nuadrupole  interactions  appear  to  be  negligible.  A value  for  the  antishielding 
parameter  X = 60  for  copper  will  be  used,  which  is  derived  from  experiments 
to  be  described  later.  This  value  is  substantially  larger  than  Sternheimer's 
value  1 + | YqqJ  = 9.7.  It  should  be  remembered  that  the  two  quantities  are 
not  strictly  comparable,  and  X will  be  larger  il  there  is  considerable  ion  core 
overlap.  In  addition  there  mav  be  a contribution  to  the  electric  field  gradient 
in  the  metal  from  distorted  valence  electron  wave  functions. 

An  upper  limit  for  the  density  of  dislocations;  which  are  assumed  to  be 
randomly  distributed,  can  be  obtained  from  the  fact  that  the  first-order 

1/2 

oroaaemng  is  less  uian  cue  aipcuar  oruaaerung  < j.  ^ , ( 


quadrupole  broadening  is  less  than  the  dipolar  broadening  ^ dipolar] 

(2.5  x 10^)  cpa.  Equation  (13)  then  yields  -an  hpper  limit  for  the  density  of 
dislc-aiio=3  in  annealed  copper  c <Clo"cm“^,  This  figure  may  be  somewhat 
too  small,  because  X's  estimate  is  rather  high.  It  is  in  fair  agreement  with 
data  derived  from  the  width  of  X-ray  diffraction  spots  on  the  assumption  of 
a random  arrangement  of  dislocations. 

Instead  of  a random  distribution  of  dislocations,  one  assume 
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dislocation  walls  and  mosaic  boundaries.  Usuallv  the  stress  field  will  be 


very  small  at  some  distance  from  the  wall  and  such  arrangements  are  not 
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detectable  by  X-ray  techniques.  It  is  possible  that  the  first  stages  of  poly- 


gonization  will  show  as  an  increase  in  g(vmax)  of  tb.e  nuclear  resonance, 


before  these  are  detectable  by  X-ray  diffraction. 


There  vras  no  evidence  lor  second-order  quadrupole  effect  on  the 
central  corqponent  in  the  cold-worked  copper.  This  is  not  unreasonable,  for 
the  calculation  of  the  broadening  of  the  central  component,  which  vanishes 
in  first  approximation  according  to  equation  (13),  in  the  stress  field  around 
a dislocation  yields  a result  which  depends  very  sensitively  on  the  cut-off 


radius  r. 

1 


The  see end  crier  interaction  decreases  sipidly  with  the  distance,, and  the 
second  moment  gives  no  satisfactory  physical  description.  The  nuclei 
are  either  so  cio&e  to  the  dislocation  that  their  contribution  to  the  central 
component  is  smeared  out  over  a large  frequency  range  and  is  thus  un- 
observable, or  else  they  are  far  enough  away  - and  this  is  the  majority  - 
that  they  produce  a nearly  unbroadened  central  line.  The  central  component 
resonance  of  nuclear  spins  a distance  r from  the  dislocation  line  will  be' 
spread  out  over  a frequency  range 


c,  extr  ~ 


21  + 3 


e4Q2  x 2 

— A. 


41  (21-1)  h vQa  r 


If  the  density  of  dislocations  is  c,  a fraction  cr  of  the  nuclei  will  be  inside 
a radius  r from  a dislocation  line. 


Fi'om  the  disappearance  of  the  satellites  in  cold- worked  copper  we 
2 4 2 

may  conclude  that  (6v)s  is  larger  than  (2.5  x 10  ) . Otherwise  th**y  would 
have  been  observable .•¥» ith the  data  of  Table  I and 
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Table  I 


1 

v cps 

for  R.^  s 55 CO  gausi 

W 

Goldschmidt 
Atomic  Radius 

■■  * 

Dipolar  Sec- 
Moment  (bvy 

ai  s 
**  i 

6.1xl06 

0 .156x10 “24 2 
cm 

1.42x10  ®cm 

4.041x10  *cm 

(3.4xl03)22 

sec 

„ .65  3 
1 

6.7xl06 

-0.145x10“  24 


-8 

1.28x10  15 

3. 60xl0~8 

(2.5xl03)2 

8 -2 

X=  60,  equation  (13)  yields  c^3  x 10  cm  for  cold-worked  copper.  On  the  other 
hand,  the  central  component  appears  to  be  unaffected.  The  central  component  is 
spread  over  a frequency  range  Av^  ®<lual  to  the  root  of  the  dipolar  second 
moment  at  a radius  r » 0.  64  10-t>  cm  according  to  equation  (14).  Assuming  that 
less  than  10  per  cent  of  the  nuclear  spins  are  closer  to  a dislocation  line,  an 

1 1 7 

upper  limit  for  the  density  of  dislocations  is  found  c <2x)0~~cm  . A value 

1 1 2 

c » 10  cm  is  compatible  with  other  experimental  evidence  for  cold-worked 
copper,  e.g. , from  stored  energy  and  X-ray  data  [ 27]  , 

In  pure  aluminum  no  decrease  in  intensity  after  cold  work  was  observed. 

TM*  observation  is  also  in  agreement  with  X-ray  data  which  shew  very. little 
broadening  of  diffraction  spots  after  cold  work  [ 28] , as  compared  to  copper. 
Presumably  the  dislocations  in  pure  aluminum  are  stacked  in  an  ordered 
arrangement  such  that  there  is  destructive  interference  of  the  stresses  at 
large  distances  from  the  dislocation  walls.  In  fact,  theory  predicts  that  the 
stress  will  drop  to  aero  according  to  the  exponential  exp(-x/A),  where  A is 
the  average  spacing  of  dislocation  lines  in  the  wall.  This  assumption  must 
be  preferred  over  the  alternative  of  self-annealing  and  over -all  lew  density 
of  dislocations,  as  the  macroscopic  mechanical  properties  of  aluminum  do 
change  after  cold  work  at  room  temperature . It  would  be  of  interest  to 
investigate  the  quadrupole  broadening  in  filings  worked  at  liquid  air  tempera- 
ture. 

The  decrease  in  the  maximum  of  the  central  component  of  the  copper 
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resonance  by  second-order  quadrupele  interaction  is  counteracted  by  a slight 
decrease  of  the  second  moment  of  the  dipolar  interaction,  when  the  satellites 
are  displaced.  This  may  seern  strange  at  first.  In  the  perfect  cubic  lattice 


of  well-annealed  copper  Van  Vleck’s  formula  equation  (1)  is  valid.  For  the 

65  65 

Cu  , resonance,  e.g.  , the  first  term  extends  over  all  other  Cu  nuclei,  the 

63  65 

second  term  over  Cu  nuclei.  When  the  four  spin  levels  of  the  Cu  nucleus 


cease  to  be  equidistant,  the  first  term  is  no  longer  correct.  In  itB  derivation 
that  part  of  the  dipolar  interaction  was  retained  in  which  a pair  of  spins  i and 


j does  not  change  its  total  rn-value,  or  terms  for  which  Am.  = Am.  = 0 or 

J 

Am.  =-  Am.  = + 1.  In  the  cold-worked  copper  part  of  the  dipolar  interact! 

1 J 3 1 1 1 

corresponding  to  transitions  like  Am.  = y—^y>  Am.  = - y have  to  be  ex- 
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interaction.  The  situation  in  which  this  interaction  vanishes  accidentally  be- 
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serious,  the  larger  the  quadrupole  interaction.  Furthermore,  only  transitions 

+ 1 


produced  by  the  external  oscillating  field  between  the  levels  - — are  of 
interest.  In  equation  (1)  transitions  between  all  levels  were  taken  into  ac- 
count. Kambe  and  Ollon  [29]  have  derived  a formula  for  the  dipolar  broaden- 


ing of  the  central  component  alone.  For  I = -j  the  proper  numerical  factor  iB 
now  instead  of  ^ in  the  first  term  and  for  I = 4 the  new  factor  10  * “ 5 
The  factor  of  the  second  term  remains  unchanged. 


Tff5X  T 


There  are  also  some  changes  in  the  establishment  of  thermal  equili- 
brium, if  the  levels  are  unequally  spaced.  In.  metals  the  important  relaxation 
mechanism  consists  of  magnetic  dipole  transitions  caused  by  interaction  with 
the  conduction  electrons.  The  set  of  differential  equations  for  the  population 
in  each  of  the  sublevels  is  solved  in  the  appendix.  The  result  is  that  in  the 
steady-state  method  of  saturation  the  maximum  of  the  imaginary  part  of  the 
susceptibility  drops  to  half  its  value  for  ^ y T^  ?(vmax)  = i.e.,the 

same  result  as  in  the  absence  of  quadrupole  broadening;  T ^ is  the  relaxation 
time  for  the  caBe  of  equidistant  levels.  The  recovery  to  thermal  equilibrium 


is,  however,  appreciably  different,  if  the  m-levels  are  unequally  spaced. 
Suppose  that  the  population  in  the  and  ---^  levels  has  been  equalized 


by  complete  saturation.  The  recovery  of  the  Boltzmann  distribution  f o .=  these 


two  levels  is  described  by  a linear  combination  of  two  exponentials  with 


characterist?  ~ times  of  T,  and  T,  The  recovery  of  the  signal  in  small 
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radiofrequency  fields  after  saturation  of  the  central  component  with  I = ■» 
is  given  by 

• / *T  fa  /rr 

“1/  A j "UW  a j 

gc(t>  = (-  4 £ -4  e + Dgc 

The  initial  recovery  rate  is  four  times  as  fast  as  in  the  absence  of  quadrupole 

5 

perturbation;  and  9 times  as  fast  for  I = ■j  . This  rate  is  observed,  e.g.,  in 
pulse  or  spin-echo  method  if  only  the  central  component  is  flipped  over. 
Considerable  care  in  the  analysis  of  the  recovery  of  the  spin  system  is  re- 
quired in  each  individual  case  . 


B.  Vacancies  and  interstitials 


The  calculations  of  the  nuclear  quadrupole  interaction  around  a static 
vacancy  or  interstitial  will  be  along  the  same  lines  as  the  calculations  for 
chemical  impurities  described  in  the  next  paragraph.  There  is  a good  reason 
not  to  treat  it  here,  as  these  effects  become  important  only  for  concentrations 
higher  than  0.1  per  cent.  Svch  high  concentrations  of  vacancies  normally 
occur  only  at  high  temperatures  close  to  the  melting  point,  where  there  is  at 
the  aame  time  a rapid  rate  of  diffusion.  In  the  case  of  severe  radiation 
damage  produced  at  low  temperature  a static  distribution  of  bigh  vacancy  and 
interstitial  concentration  could  be  obtained.  For  the  treatment  of  this  case  we 
refer  to  the  next  section.  Here  the  time-dependent  quadrupole  effect  in  a 
crystal  with  an  appreciable  vacancy  concentration  at  elevated  temperatures 
is  discussed.  We  can  describe  the  diffusion  as  a random  motion  of  the 
vacancies  through  the  lattice.  Assume  for  the  sake  of  simplicity  that  the 
electric  field  is  q whenever  a nucleus  is  adjacent  to  a vacant  site  and  zero 
elsewhere.  A more  rigorous  treatment  would  obviously  include  a summation 
over  a distribution  of  q-values.  The  probability  for  a nucleus  to  be  next  to 
a vacancy  is 


-E 


p(0)  = - ze 


vac 


ykT 


where  z is  the  number  of  neighboring  sites  and  Evac  is  the  energy  required 
for  the  formation  of  a Schottky  vacancy.  The  probability  that  a given  nucleus 
is  still  next  to  the  same  vacwicy  after  a time  t is  exp(-t/T  ) where 


i..iiifc)»wiaiatw«ww» 


TR199 


-23- 


r'1  - (I  - l)vat  e 
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is  the  jumping  Ire qaency  of  the  vacancy  to  any  of  the  z-1  other  neighboring 
positions.  The  correlation  function  of  the  time-dependent  electric  Held 
gradient  is  consequently 


2 

= q ze 


-E  /kT  -t  '/T 
vac 


and  the  quadripole  contribution  to  the  transition  probability  or  the  inverse 
nuclear  relaxation  time  is  consequently  proportional  to 
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In  this  most  important  case  of  fast  diffusion  the  relaxation  by  quadrupole 
coupling  contains  a different  combination  of  energies  in  the  exponent  than  the 
dipolar  interaction,  which  pontains  {Ears^  + Ear,t)  as  does  the  diffusion  constant. 
If  the  quadrupole  interaction  is  dominant  at  high  temperatures,  E and  E . 
can  be  determined  separately  with  resonance  methods.  The  condition 
/T<^Cvo  WlO  sec  will  usually  be  satisfied  only  quite  close  to  the  melting 
point.  At  lower  temperatures  where  vq~  * , the  quadrupole  interaction  will 
usually  be  negligible  compared  to  magnetic  interactions.  The  same  treatment 
holds  for  interstitials  except  that  z/z-1  is  then  replaced  by  unity  in  equation 
(15).  No  experimental  evidence  of  the  effect  described  by  this  formula  is 
known.  It  could  occur  both  in  metals  and  insulators. 

The  quadrupole  interaction  in  liquids  and  its  contribution  tc  the  relaxation 
mechanism  is  well  established  [ 30,31]  . One  could  consider  this  as  a limiting 
case  of  the  foregoing  in  which  the  concentration  of  vacancies  is  large  and  in- 
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dependent  of  file  temperature.  depends  on  the  temperature  only  through 
the  jumping  time  T which  is  proportional  to  (T|  / T)  where  Y|  is  the  viscosity. 

C.  Foreign  atoms 

The  influence  of  foreign  atoms  on  the  electric  field  gradient  is  two- 
fold. une  can  distinguish  between  charge  and  size  effects.  In  the  first  place 
one  can  imagine  an  ion  in  an  ionic  crystal  replaced  by  another  of  the  same 
size  but  differing  by  an  electronic  charge  e.  A vacancy  will  also  act  as  an 
extra  charge,  equal  and  opposite  to  the  charge  of  the  ion  normally  occupying 
the  lattice  site.  Some  distance  r away  from  Mie  impurity  atom  the  quadru- 
pole  interaction  with  a nucleus  will  cause  a change  in  the  resonant  frequency 
between  the  m and  m - 1 spin  level 


h Av  = (2m- 1) 


3 e2Q 
■xvi  acv  "i  r 
oi\ 


1 +l^aJ  2 ,, 

y (3  COS  0 -1) 


where  c is  the  dielectric  constant,  9 tne  angle  between  the  radius  vector  r 
and  the  external  magnetic  field,  while  (1  + IYqJ)  the  antisfcielding  factor 
calculated  by  Foley  and  collaborators.  In  a metal  an  impurity  atom  of  the 
same  atomic  volume  but  different  valency  will,  in  the  free  electron  approxi- 
mation, produce  a screened  Coulomb  field  [12]  around  it  with  screening  para- 
meter K,  and  the  change  in  splitting  of  the  spin  levels  due  to  the  quadrupole 
interaction  will  be  give*  by 


hAv  = -ii*"1: U.2*  Q.  (1  + Kr  +yKZr2)  e'Kr  (1  <-|  y |)(3coB2e-i) 
81(21  -ljr5  * 00 

Take  the  square  of  the  displacement  and  average  over  a sphere  of  volume 
■jirr5  - C~\  where  C is  the  impurity  concentration,  to  obtain  the  second 
moment  of  the  broadened  satellites.  The  first  moment  is  zero  because  of 
the  angular  dependence.  In  an  ionic  crystal  the  second  moment  will  be  given 
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This  result  depends  very  sensitively  on  an  inner  cut-off  radius  r. . The 

-9  1 

fourth  moment  would  vary  as  r^.  There  are  a few  very  large  contributions 

to  the  moments  from  nuclei  in  the  immediate  neighborhood  of  the  impurity, 

while  the  majority  causes  only  a relatively  small  broadening.  Under  these 

circumstances  a useful  physical  description  is  again  obtained  by  considering 

the  nuclei  in  a spherical  shell  of  radius  r around  the  impurity.  The  signal 

from  these  nuclear  spins  is  spread  out  over  a frequency  interval  given  by 

equation  (16),  if  v/e  substitute  3 for  the  range  of  values  of  3 cos^0  - 1.  When 

this  frequency  range  exceeds  a certain  critical  value,  which  depends  on 

the  experimental  sensitivity,  we  may  say  that  the  contribution  of  these  nuclei 

is  unobservable,  because  it  is  spread  '-rnt  over  such  a broad  frequency  range. 

A practical  value  of  thi  s limit  i s determined  by  the  signal -to -noise  ratio  and 

is  about  25  kc/s  for  the  experiments  on  the  copper  and  aluminum  resonance 

described  in  this  paper.  A critical  radius  r is  defined  bv  taking  Av  in 

cr 

equation  (16)  of  the  same  order  of  magnitude  as  the  dipolar  broadening.  The 
nuclei  inside  this  sphere  contribute  atytofhe  unobservable  wings  of  the  reso- 
nance signal;  nuclei  outside  the  sphere  give  an  essentially  unbroadened  signal. 
Clearly  this  reduction  of  the  quadrupole  interaction  to  an  "all  or  nothing" 
effect  is  a great  oversimplification.  Its  only  merit  is  its  usefulness  in 
discussing  existing  experimental  data  on  the  maximum  absorption  as  a func- 
tion of  impurity  content.  When  more  precise  experimental  data  on  the 
complete  line  shape  become  available,  the  analysis  can  be  made  more 
complete.  The  contribution  of  the  nuclear  spins  in  each  spherical  shell 
has  the  shape  shown  in  Fig.  9.  These  contributions  should  be  added  to  ob- 
tain the  resultant  line  form. 


In  addition  to  the  electric  field  gradient  produced  by  the  extra  charge 
of  vacancy,  interstitial  or  foreign  atom,  there  is  the  effect  from  the  stress 
field  around  the  imperfection.  The  strain  deformation  around  a spheilcal 
inclusion  in  an  elastic  medium  falls  off  as  the  inverse  cube  of  the  distance 
[ 32] . If  the  foreign  atom  and  solvent  atom  have  radii  and  aQ  respectively, 
the  relative  displacement  of  two  ions  a distance  r from  the  impurity  has  the 
order  of  magnitude 


a / , 3 -3 

Aa  = (a.  - a ) a r 
c ' i o'  o 


Introducing  the  same  antishielding  factor  X as  in  the  discussion  of  ion  core 
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distortion  around  a dislocation,  we  obtain 


18  e 


2 


hAv  = + 


Q(a. 


ao> 


“3 — 

r a 


(3  cos*  8 - 1) 


(17) 


o 

This  shift  has  to  be  added  to  the  result  (16)  before  calculating  the  second 
moment  of  the  satellites  by  the  quadrupole  interaction  in  an  ionic  crystal. 


It  would  be  of  interest  to  study  these  charge  and  size  effects  in 
detail  with  controlled  impurities  in  ionic  crystals . Reif  has  very  recently 

observed  such  quadrupole  effects  on  the  bromine  resonance  in  AgBr,  con- 

-f-f  + 

taining  Cd  ions.  Presumably  there  will  be  a corresponding  number  of  Ag 

vacancies  to  retain  over-all  neutrality.  At  low  temperatures  the  Cd++  ion 

and  the  vacancy  will  be  associated  [33j  . At  high  temperatures  they  are 

xiiucpendcnt  impurities.  The  quadrupole  broadening  in  such  a crystal  will 

undoubtedly  change  with  temperature  and  the  dissociation  process  could  in 

principle  be  followed  in  some  detail.  While  this  suggestion  was  written  in 

the  manuscript,  the  author  was  unaware  that  Reif  had  already  discovered 

such  effects  and  their  explanation  had  already  been  given  by  Cohen  and 

Reif  [34]  . 


C.  1 First-order  quadrupole  effects  in  metals 

8 - 1 

In  a metal  the  electric  field  is  so  strongly  shielded  (K  = 2.10  cm  ) 
that  the  ion  deformation  by  "trains  gives  the  only  contribution  to  quadrupole 
broadening.  Equation  (17)  should  be  used.  In  Fig.  14  the  maximum  ab- 
sorption  of  the  satellites  of  A -U  in  an  A -6  Mg  alloy  observed  by  Rowland 
[ 35]  is  plotted  as  a function  of  Mg -concentration.  Cold-worked  filings  dis- 
play at  first  the  combined  effect  of  impurities  and  dislocations.  The  fact 
that  in  cold-worked  impure  aluminum  the  dislocations  are  not  annealed  out 
at  room  temperature  is  in  agreement  with  other  metallurgical  observations. 
According  to  Cottrell  [36]  the  dislocations  are  locked  to  an  "atmosphere" 
of  impurity  atoms.  On  annealing  at  higher  temperatures  the  concentration 
of  dislocations  becomes  too  small  to  give  an  observable  contribution.  Tne 
first-order  quadrupole  effect  of  the  impurities  is  then  retained.  We  may 
subtract  from  the  total  observed  absorption  g(*max)  the  central  component 

contribution  g (v  ) and  plot  the  maximum  intensity  of  the  satellite s 
6c  max'  r ' 
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Fig.  13.  The  line  shape  of  the  At  resonance  in  aluminum  alloy  with 
0.64  atomic  per  cent  magnesium,  (a)  after  filing  at  room  temperature, 
(b)  after  filing  and  87  hour":  at  250°C,  (c)  after  filing  and  2 hours  at 480°C, 
(d)  the  absorption  in  pure  aluminum,  included  for  comparison.  The  peak 
intensity  increases,  as  the  density  of  dislocations  decreases  with  anneal- 
ing. The  remaining  reduction  of  (c)  as  compared  to  (d)  is  the  effect  of 
the  Mg-impurities  alone.  (Rowland). 


a*-  *t.war  .•«—*». 


C,  ATOMIC  FRACTION  Mg 

Fig.  14.  The  peak  intensity  of  the  resonance  in  A(~Mg  in 

annealed  alloy  = as  a function  of  the  Mg -concentration  at  room  tem- 
perature. The  satellite:  peak  intensity  decreases  rapidly.  At  4 
atomic  per  cent  Mg  only  a contribution  of  the  central  component 
remains.  This  component  extrapolates  to  0.257  in  pur e aluminum, 
the  theoretical  value  for  I ~ 5/2.  The  alloys  with  high  Mg-concen- 
tration  were  obtained  by  quenching  the  powdered  specimen  from 
high  temperatures,  where  the  solid  solubility  is  high.  (Rowland). 


r ■■  '■ 


?.Hv>  iU  L&  H * •.rVii 


TR199  -27- 

g (v  ) in  the  annealed  alloy  as  a function  of  the  Mg -concentration  on  a 
sa-  max'  7 6 

double  logarithmic  scale.  The  result  in  a straight  liue.  This  finds  a 
natural  explanation  in  terms  of  the  "all  or  nothing"  language  introduced 
earlier.  If  there  is  no  Mg -atom  inside  the  critical  radius  the  aluminum 
nucleus  will  contribute  to  the  satellite  signal;  otherwise  not.  The  Mg-atoms 
are  supposed  to  be  distributed  at  random  over  the  lattice,  which  is  certainly 
correct  at  these  low  concentrations.  Sach  lattice  has  a probability  C to  be 
occupied  by  an  Mg  atom,  where  C is  the  relative  Mg  concentration.  The 
number  of  lattice  sites  inside  the  critical  radius  is 

4t  3 ; 

N = r /v_ 

3 era 

The  probability  that  an  A l nucleus  has  no  Mg  atom  inside  the  critical  distance 

N 

r is  consequently  (1-C)  and  the  maximum  intensity  of  the  satellite  signal 
cr  c 

for  I = -j,  normalised  per  A A atom 


a ( v ) = g 
ss  max7 


s.O^maxX1-^ 


Experimentally  it  is  found  that  N = 130.  Since  N is  so  large  the  cut-off  is 
not  as  sharp  as  we  have  assumed.  It  is  somewhat  surprising  that  the  ex- 
perimental results  can  be  so  well  described  in  this  fashion. 

The  introduction  of  zinc  in  the  aluminum  lattice  gives  a reduction  of  the 

satellites,  shown  in  Figs.  (15)  and  {16),  which  are  well  described  by  N = 98. 

The  cut-off  radius  lies  between  the  sixth  and  seventh  near  neighbor  in  the 

face-centered  cubic  lattice.  The  ratio  of  quadrupole  interaction  at  these 

3 

two  sites  is  q^/q^  = (ry/r£)  = This  fraction  is  not  large  and  the  cut- 

off is  not  very  sharp.  If  the  experimental  precision  were  increased  - e.g.  , 
by  inc reading  the  absorption  intensity  through  cooling  to  liquid  helium 
temperature  - a more  precise  analysis  of  the  strain  field  could  be  made 
Wings  of  the  satellite  away  from  the  central  peak  have  already  been  observed, 
but  a precise  evaluation  of  the  line  shape  has  not  yet  been  made. 

It  is  not  surprising  that  the  effective  cut-off  radius  for  Zn  is  smaller, 
as  the  Zn  atom  has  more  nearly  the  same  size  as  the  A*t-atom  than  mag- 
nesium. Since  the  gradient  at  the  cut-off  radius  must  be  the  same,  it  follows 
from  equation  (17),  that 
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This  is  in  poor  agreement  with  tabulated  values  of  Hume-Rothery  [ 36) , which 
give  3.5  for  this  ratio. 

C.2  Second-order  quadrupole  effects 

The  broadening  of  the  central  component  can  be  analyzed  along  similar 
lines.  Since  this  interaction  will  drop  off  with  the  inverse  sixth  power  of  the 
distance  from  the  impurity,  the  second  moment  of  this  interaction  averaged 
over  the  whole  crystal  will  give  even  less  relevant  information  than  for  the 
first-order  effect.  The  second  moment  would  be  proportional  to  the  inverse 
ninth  power  of  inner  cut-off  radius.  The  nuclear  spins  at  a distance  r from 
the  impurity  will  give  a contribution  to  the  central  component  which  has  the 
shape  displayed  in  Fig.  10  with  a distance  between  the  extremes  given  by 


21  + 3 


4 

e Q 
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Assuming  that  the  inverse  cube  law  for  the  strains  holds  also  in  the 
immediate  neighborhood  of  the  impurity,  this  frequency  range  at  first, 
second  ana  third  nearest  neighbors  in  a face-centered  cubic  lattice  has  the 
ratios  27:  3.37:  1.  It  is  therefore  quite  possible  that  this  range  for  third 
nearest  neighbors  is  only  a fraction  of  the  dipolar  broadening,  while  it  is 
appreciably  larger  than  the  dipolar  interaction  for  nearest  and  next  nearest 
neighbors.  The  nuclear  spins  which  have  an  impurity  in  one  of  these  z 
neighboring  sites  would  then  give  an  observable  effect,  while  other  nuclei 
would  contribute  to  an^>thervjistf  unchanged  central  component.  In  other 

f 

words  a sharp  cut-off  may  now  be  expected  with  more  justification  than  in 
the  satellite  case.  The  maximum  intensity  of  the  central  component  if  I 
is  an  odd  half-iuteger  is  consequently  given  by 


g { v ) = g ( v ) ( 1 - C) 
Bc  max7  *c,o'  max'  ' ' 


where  C is  the  relative  impurity  concentration  of  the  solute  atoms.  For 
a face-certered  cubic  lattice,  z may  assume  one  of  the  values  0,  12  or  18, 
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Fig.  15.-  The  peak  intensity  of  the  A-£^  resonance  in  a£  -Zn  alloys  at 
room  temperature,  at  6.18  Mc/f;.  The  satellites  are  washed  cut,  the 

central  component  remains  unchanged.  (Rowland), 
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Fig.  16.  A double  logarithmic  plot  of  the  peak  intensity 
of  the  satellites  alone 

^vmax^  9 26 _ ^8^vmax' 

g (v  ) 3?  25  v.  (v  T 

°o  max  -»,o  max 

for  the  Ai  ^ resonance  in  A^-Zn  alloys  as  a function  of 
the  Zn-concentration.  The  straight  line  has  a slope  of 
98.  The  data  have  been  obtained  from  Fig.  15  by  sub- 
tracting the  effect  of  the  central  component.  (Rowland). 
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etc.  and  iu  a body-centered  cubic  lattice  0,  6 or  14.  The  data  on  the 

central  component  are  moat  conveniently  taken  in  cold-worked  samples. 

The  contribution  of  satellites  to  the  signal  is  then  effectively  eliminated  by 

the  large  density  of  dislocations  so  that  no  correction  for  them  is  necessary. 

In  this  way  the  data  in  Fig.  17  for  the  central  component  in  a-brass  have 

been  obtained  [ 1 1]  . Indeed  it  was  found  that  z = 18  for  the  copper  resonance 

27 

in  a-brass,  and  z = 0 for  the  A'u  ‘ resonance  in  the  aluminum  zinc  alloy. 

These  numbers  depend  on  the  resonance  frequency  vq.  In  very  high 

external  magnetic  fields  the  second-order  cruadrupole  interaction  will  tend 

pZ7 

to  zero, and  so  does  the  critical  radius.  For  lower  frequencies  the  Av 
resonance  in  the  aluminum  zinc  system  may  give  z = 12.  It  would  be  inter- 
esting to  follow  nearly  discontinuous  changes  in  z,  as  »o  is  increased  or 
decreased.  At  very  low  fields  the  quadrupple  interaction  might  extend  to 
third  nearest  neighbors  etc.  and  a detailed  study  of  this  kind  will  give  in- 
formation about  the  stress  field  in  the  immediate  neighborhood  of  the  im- 
purity and  will  show  any  deviation  from  the  continuum  theory  which  might 
exist  on  an  atomic  scale. 

From  the  first-order  and  second-order  broadening  values  for  the 
product  X(a.  - aQ/a)  may  be  derived.  A value  (a^  - aQ/a)  = 10  per  cent 
will  be  used  for  a-brass,  and  other  pertinent  data  appeared  in  Table  1. 

If  we  assume  that  the  quadrupolar  spread  in  frequency  Av  is  equal 

to  the  dipolar  broadening  ^(6v)  y for  third  nearest  neighbors,  and 
3.37  times  as  large  for  second  nearest  neighbors,  we  find  X = 60  for  the 
copper  lattice.  This  estimate  is  not  very  reliable,  as  the  elastic  contiuum 
theory  is  not  valid  for  near  neighbors.  It  is  much  larger  than  the 
Sternheimer  antishielding  factor.  This  may  be  correct  as  the  effects  from 
ion  cere  overlap  will  be  severe  in  the  close -packed  copper  lattice.  This 
rather  large  value  of  X.  has  been  used  previously  to  estimate  the  density  of 
dislocations  in  copper. 

A more  reliable  value  X in  aluminum  may  be  obtained  from  the  first- 
order  quadrupcle  effect  in  aluminum  alloys.  The  critical  radius  for  satellite 
broadening  in  the  A&-Zn  alloy  for  N = 98  is  rcr  = 6.5  x 10"°  cm.  At  this 
radius  the  mean-square  broadening  of  theinner  satellites  is  assumed  to  be 
twice  the  dipolar  second  moment.  With  (a^  - aQ)/a0  - 4 per  cent  for  zinc 
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in  aluminvm,  equation  (17)  yield#  X.  = 30:  Probably  (a^  - a0)/a  should  be 

taken  somewhat  larger  and  X correspondingly  smaller . The  fact  that  no 

second-order  quadrupole  effect  is  observed  in  this  alloy  is  consistent  with 

the  value  X(a.  - a )/a-«  1.2.  The  broadening  of  the  central  component  for 
l o 2 

nearest  neighbors  is  Av  = 2.4  x 10  cps,  or  less  than  the  dipolar 

c (Cxir 

broadening.  The  margin  is  small,  however. 

In  the  A^-Mg  alloy  the  strains  in  corresponding  lattice  points  are 
larger  by  a factor  1.3,  and  in  this  case  the  broadening  of  the  central  com- 
ponent for  nearest  neighbors  is  barely  noticeable.  The  observed  line  in  a 
quenched  alloy  containing  14  atomic  per  pent  Mg  reproduced  in  Fig.  18, 
is  a superposition  of  a line  shape  ofFi&.  10  for  A-6  nuclei  which  have  an 
Mg-atom  as  nearest  neighbor  and  an  undistorted  contribution  of  other 
At  nuclei.  The  asymmetry  is  much  more  pronounced  in  the  derivative 
curve,  which  was  actually  recorded.  Clearly  the  second-order  broadening 
from  nearest  neighbors  is  small.  If  we  extend  the  strain  theory  to  the 

neighbors,  the  second-order  quadrupole  interaction  should  be  a factor 

2 2 

^aMg  " ao^aZn  ~ ao>  = larger  *or  magnesium  in  aluminum  than  for 

zinc  in  aluminum  or 

Av  ^ = (1.3)^  x 2.4  x 10^  = 4.1  x lO^cps 

This  is  somewhat  larger  than  the  dipolar  width  and  a small  effect  is  ob- 
served in  the  magnesium  alloy,  whereas  no  effect  was  observed  in  the 

zinc  allov  at  the  freouencv  of  observation  v = 6.1  x 10^  cos.  Thus  the 
• * - o 

magnitude  of  the  second-order  interaction  is  determined  rather  precisely. 

No  noticeable  shift  of  the  resonance  maximum  was  found  in  the 
aluminum  and  copper  alloys.  This  is  in  agreement  with  the  discussion  of 
the  Knight  shift  in  alloys  illustrated  in  Fig.  4.  Only  those  copper  and 
aluminum  nuclei  which  are  a few  interactomic  distances  away  from  the  solute 
atoms,  contribute  tc  the  maximum  of  the  resonance  absorption.  At  these 
positions  the  product  of  density  of  states  and  density  of  the  one -electron  wave 
function  remains  unchanged.  The  resonance  of  nuclei  adjacent  io  he  impurity, 
which  could  .show  a change  of  the  Knight  shift,  is  unobservable  due  to  quadru- 
pole interaction.. 
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18.  The  dashed  line  is  the  intensity  which  could  be  expected 
in  an  ordered  arrangement  where  ail  zinc  atoms  are  contained 
in  one  simple  cubic  sublattice. 
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C.  3 Precipitation 

The  values  of  \ derived  for  the  copper  and  aluminium  lattice  were 
used  previously  to  evaluate  the  density  of  dislocation  lines  in  these  cold- 
worked  metals.  It  has  been  shown  how  annealing  of  metals  and  alloys  is 
generally  accompanied  by  an  increase  in  the  maximum  magnetic  resonance 
absorption,  as  the  density  of  dislocations  decreases.  This  effect  of  tempera- 
ture treatment  should  be  sharply  distinguished  from  the  precipitation  of  new 
ph*s^».  The  latter  will  frequently  lead  also  to  an  increase  of  the  resonance 
signal. 

If  a precipitate  is  formed  the  remaining  matrix  will  contain  less  solute 
atoms  and  the  signal  g(v  max)  ought  to  increase  rapiily,  ae  the  impurity  con- 
centration C in  the  original  matrix  decreases.  At  firaf  this  effect  may  be 
obscured  by  the  strains  set  up  by  the  precipitation  process.  The  intensity 
“*Vmax'  ~e  ““  a function  of  thermal  treatment  and  interpreted 

in  terms  of  these  two  effects.  No  detailed  analyses  have  yet  been  made, 
but  the  occurrence  of  precipitation  in  a quenched  14  per  cent  Mg  solution 
in  A-£  after  heating  for  five  days  at  250°C  is  shown  in  Fig.  18.  The  asymme- 
try due  to  secondvorder  perturbation  of  the  central  component  disappears, 
ac  at  the  end  of  the  heat  treatment  the  aluminum  matrix  contains  only  3 per 
cent  or  4 per  cent  Mg  atoms  and  most  aluminum  nuclei  do  not  have  an  Mg 
atom  as  neighbor.  If  the  quadrupole  interaction  had  been  larger  (v  smaller) 
the  effect  would  have  been  more  pronounced.  The  increase  in  intensity  after 
precipitation  and  annealing  should  be  especially  striking  if  the  final  matrix 
would  contain  so  few  solute  atoms  that  the  satellites  would  reappear. 

C.4  Order-Disorder 

■■■  — ■ .i  i ■■  * 

Since  the  central  component  signal  depends  so  sensitively  on  the 
atomic  configuration  in  the  immediate  neighborhood  of  the  nucleus  with  the 
quadrupole  moment,  8c(vmax)  maY  he  used  to  investigate  short-range  order 
phenomena.  Equation  (20)  derived  under  the  assumption  of  completely 
random  solution  has  to  be  modified,  if  short-range  order  is  present.  We 
shall  first  discuss  the  case  of  the  maximum  short-range  order  which  could 
occur  in  the  face-centered  cubic  lattice  of  a-brass.  This  is  the  situation 
at  absolute  zero  of  temperature,  if  the  reaction  times  would  be  sufficiently 
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rapid  for  it  to  be  established.  Clearly  this  situation  will  never  be  attained 
completely,  but  this  theoretical  case  will  then  show  what  deviations  from 
equation  (20)  ran  at  most  be  expected.  Only  interaction  between  nearest 
neighbors  is  assumed  and  V + V -2V  >0  in  a brass.  Pre- 

cipitatlon  of  zinc  would  occur,  if  the  opposite  sign  were  used.  The  Zn 
atoms  will  tend  to  st ey  epert  sufficiontiy  low  temos r sturc s the 

equilibrium  situation  will  be  one  in  which  all  Zn’s  belong  to  the  same 
simple  cubic  sublattice  of  the  fee  lattice;  in  this  sublaiiice  the  distribution 


of  Zn-atoms  will  be  random.  The  model  could  be  refined  to  take  into  account 
interaction  between  next  nearest  neighbors  which  are  neighbors  in  the  same 
oublattlce.  A general  theory  which  covers  this  case  will  be  mentioned  shortly. 
If  the  relative  concentration  of  Zn  is  C<0.25,  the  relative  concentration 
on  one  subiattice  is  4G.  The  other  three  sublattices  are  pure  copper.  Cal- 
culation of  the  probability  that  a copper  atom  has  cubic  symmetry  for  the 
configuration  of  its  18  neighbors  gives  for  the  maximum  of  the  shape  function 
per  copper  nucleus  in  this  model 


g (y  )/ a ( v ) 

6c'  max'  ®c, o'  max' 


3C 

*TT^o 


(1  - 4C)  + 


1 -4C 

XT-rcy 


(i-4c)°  * (4C y 


(21) 


The  first  term  is  the  contribution  to  the  intensity  of  the  central  component 
from  the  three  pure  copper  sublattices.  The  last  term  is  from  the  remaining 
copper  atoms.  There  always  have  other  copper  atoms  as  nearest  neighbors; 
they  contribute  to  the  signal  if  the  six  next  nearest  neighbors  are  either  all 
copper  or  all  zinc  atoms.  The  result  is  indicated  as  the  dotted  line  in  Fig.  17 
For  low  Zn  concentrations  the  result  is  the  same  as  for  the  random  solution. 
At  low*concentrations  the  Zn  atoms  are  far  away  from  each  other,  regardless 
of  whether  the  atate  is  ordered  or  disordered.  Each  Zn  atom  excludes  18 
copper  atoms  from  giving  a significant  contribution  and  the  intensity  drops 
as  1-18C  in  this  region.  At  higher  concentrations  there  are  fewer  Cu  nuclei 
contributing  to  the  signal.  Since  the  Zn  atoms  tend  to  stay  apart,  there  are 
more  copper  nuclei  with  a Zn  atom  in  their  neighborhood  than  in  the  random 
situation.  The  peculiar  secondary  maximum  near  20  per  cent  concentration 
is  caused  by  configurations  with  six  Zn-atoms  as  next  nearest  neighbors. 

In  principle  careful  measurements  of  the  intensity  of  the  magnetic  resonance 


i 
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Fig.  18.  (a)  The  shape  of  the  A*27  resonance  in  a quenched  14  per  cent 
Mg  alloy.  The  line  is  slightly  asymmetric  by  incipient  second-order 
quadrupole  perturbation,  (b)  the  shape  in  the  same  sample  after  heat 
treatment.  A Mg-rich  phase  precipitates,  the  signal  of  the  A-fc -nuclei 
in  the  remaining  matrix  increases.  The  Mg  content  is  still  largerthan 
4 per  cent,  as  the  contribution  from  the  satellites  remains  negligible. 
(Ro  viand). 
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of  copper  in  a-brass  with  Zn  concentrations  between  10  and  20  per  cent 
would  give  definite  information  about  the  occurrence  of  ordering.  Such 
measurements  could  be  carried  out  with  sufficient  accuracy  at  liquid  helium 
temperatures.  For  the  completely  ordered  state  at  25  per  cent  Zn  the  in- 
tensity drops  to  zero. 

Although  this  situation  cannot  be  realized  with  a-brass,  the  same 

geometry  is  attained  in  Cu,  Au.  The  Cu  atoms  in  this  ordered  cubic 

alloy  are  not  center*  of  cubic  symmetry.  They  have  8 Cu  atoms  and  4 Au 

atoms  as  nearest  neighbors.  Therefore  the  copper  resonance  will  be  un- 

obeervable  in  this  ordered  alloy.  For  the  investigation  of  hhort-range  order 

197 

near  the  composition  Cu^Au,  the  Au  resonance  should  be  used.  The  Au 
atoms  have  a cubic  environment  in  che  ordered  state.  For  lower  Au  con- 
centrations (between  0-15  per  cent)  the  copper  resonance  will  be  a suit- 
able tool  to  investigate  the  short-range  ordering  in  the  copper-gold  system. 

The  theory  of  short-range  order  best  adapted  to  our  problem  is  that 
of  Cowley  [ 38]  , who  introduces  a short-range  order  parameter  for  every 
lattice  site.*  Take  the  origin  at  a copper  nucleus.  The  probability  that  a 
site  with  coordinates  ■{,  m n relative  to  this  origin  is  also  a copper  nucleus  is 


where  a^mn  *8  the  order  parameter  as  defined  by  Cowley  and  C is  the 
relative  Zn  concentration.  Cowley  shows  how  to  obtain  a set  of  equations 
for  cu,  the  value  of  in  the  i^  shell  around  the  origin,  for  any 

temperature  and  concentration.  He  gives  explicit  results  for  a^,  a^  and 
for  Cu^Au  as  a function  of  temperature.  These  agree  well  with  experi- 
mental X-ray  data.  The  a£mn  are  the  coefficients  in  the  three-dimensional 
Fourier  series  expressing  short-range  order  "scattering  power”  as  a.  function 
of  reciprocal -lattice  coordinates.  In  terms  of  Cowley's  parameters  and 
the  intensity  8c(vmax)  copper  resonance  in  a-brass  becomes 


'Ac, 


a (v  )/g  (v  ) = (1  -C+Ca.)12  (1  -C  f Ca,)6 
®cx  max'/  Bc,o'  max'  ' 1#  ' 2' 


o'  max' 


*The  author  is  indebted  to  Professor  H.  Brooks  for  calling  his  attention  to 
Cowley's  theory. 
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Here  the  possibility  that  all  twelve  nearest  and/or  six  next  nearest  neighbors 
are  zinc  has  been  neglected.  This  method  can  be  extended  in  a straightforward 
fashion  to  cases  where  other  numbers  of  neighbors  are  involved. 


C.  5 Non-  cubic  crystals 

Finally  a few  remarks  must  be  made  about  the  influence  of  dislocations 
and  other  imperfections  on  the  resonance  in  .nan-cubic  single  crystals  and  on 
pure  quadrupole  transitions. 


An  insulating  ncn~cubic  single  crystal  has  a definite  advantage  over 
cubic  samples  in  the  study  of  quadrupole  effects,  because  the  central  compo- 
nents and  the  satellites  are  well  separated  in  frequency  even  in  the  perfect 
crystal.  The  broadening  by  imperfections  can  be  analyzed  for  each  line 
separately.  The  results  obtained  by  Pound  [19]  on  a single  crystal  of 

indicate  that  the  satellites  are  indeed  broader  than  the  central  com- 
ponent, and  the  outer  pair  is  again  broader  than  the  inner  pair  of  satellites. 
The  presence  of  dislocations  or  impurities  presumably  causes  a distribution 
of  q-value  s around  the  single  value  for  the  perfect  crystal  and  the  influence  on 
each  component  line  can  be  analyzed  in  the  manner  indicated  above.  Using 
equation  (13)  and  taking  for  X in  A-fc.,0,  the  same  value  as  in  aluminum,  the 

density  cf  dislocations,  if  these  arc  assumed  to  be  the  cause,  is  csiimai 

7 -2 

from  Found’s  recordings  to  be  c^-S.lO  cm 


It  is  not  certain,  however,  that  the  broadening  is  only  caused  by 
a random  distribution  of  dislocations  in  this  non-cubic  single  crystal,  as 
opposed  to  a previous  conclusion  valid  for  powder  samples.  Slightly  different 
orientations  of  the  axis  of  symmetry  in  various  subgrains  of  the  crystal  due 
to  dislocation  walls  will  also  produce  a broadening.  The  magnitude  of  the 
spread  in  resonance  frequencies  caused  by  a distribution  A0  of  the  angle  6 
i'4  obtained  by  differentiation  of  equation  (12).  A spread  A0  of  0.5  degree 
would  give  the  required  magnitude  of  tic  broadening,  which  should  vary  as 
sin  2 9.  This  angular  dependence  is  not  observed  experimentally,  and  there- 
fore random  dislocations  are  more  likely  to  be  the  broadening  agent. 


In  zero-field  pure  quadrupole  spectra  the  presence  of  imperfections 
.will  cause  a first-order  broadening  of  the  lines.  The  order  of  magnitude 
will  be  the  same  as  that  given  by  equations  (13)  and  (16)  respectively. 


i 
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Found  [2]  has  observed  such  broadening  due  to  cold  work  lor  the  i 

5 3 

quadrupole  resonance  in  a.  molecular  crystal  of  iodine.  The  m,  = + •»-*•+•*• 

3 j 1 c c 

transition  is  broadened  twice  as  much  as  the  = + 2~^+  *2  transition. 


D.  Electrons 

The  quadrupole  interaction  of  a nucleus  with  the  conduction  electrons 
does  not  give  rise  to  prominent  effects.  In  a non-cubic  crystal  the  asym- 
metric charge  distribution  of  all  conduction  electrons  is  incorporated  in  the 
distribution  of  all  charges  in  the  crystal  to  give  the  electric  field  gradient 
at  the  nucleus,  and  it  cannot  be  observed  separately.  The  asymmetric 
distribution  of  the  conduction  electrons  near  the  Fermi  level  gives  rise 
to  an  anisotropic  Knight  shift.  The  qp  of  only  theee  electrons  may  be  ob- 
served if  a nucleus  with  a spin  I = ^ is  present. 

The  quadrupole  interaction  will  have  a small  effect  on  the  relaxation 

time  even  in  cubic  crystals.  An  electron  with  a wave  vector  t“  G near  the 

Fermi  level  can  flip  the  nuclear  spin  not  only  by  the  magnetic  dipolar 

interaction  discussed  previously  but  also  by  quadrupole  interaction,  since 

it  has  an  asymmetry  in  its  wave  function  and  individual  electrons  will 

contribute  independently  to  the  relaxation  mechanism.  As  the  hyperfine 

interaction  will  always  have  a predominantly  magnetic  character  the  reduction 

of  T.  by  quadrupole  interaction  will  be  small.  Nevertheless,  there  is  some 
1 19  ?! 

evidence  from  the  relaxation  time  [16]  of  the  isotopes  Ga  and  Ga  in 

8 5 87 

metallic  gallium  and  Rb  and  Rb  in  rubidium  that  quadrupole  interaction 
contributes.  The  relaxation  times  for  these  isotopic  pairs  are  net  inversely 
proportional  to  the  square  of  their  gyromagnetic  ratios.  The  relaxation  time 
is  relatively  shorter  for  the  isotopes  with  the  larger  quadrupole  moment. 


E.  Phonons 

The  lattice  vibrations  will  produce  a time-varying  gradient  of  the 

electric  field  even  in  cubic  crystals.  The  thermal  relaxation  by  quadrupole 

interaction  in  solids  has  been  demonstrated  in  an  unambiguous  fashion  by 

Pound  [19]  in  a single  crystal  of  sodium  nitrate.  His  publication  should  be 

consulted  for  the  elegant  method  of  double  resonance  with  simultaneous 

transitions  at  two  different  frequencies  between  the  four  unequally  spaced 
23 

levels  of  the  Na  spin.  In  many  ionic  crystals  t.he  quadrupole  relaxation 
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is  more  important  than  the  magnetic  interaction  via  paramagnetic  im- 
perfections. The  latter  is  the  only  mechanism  to  give  reasonable  re- 
laxation time  for  I = -j.  In  molecular  crystals  large  quadrupole  inter- 
actions may  occur  which  may  also  determine  the  relaxation  time.  The 
quadrupole  interaction  in  the  presence  of  molecular  rotation  or  tunnelling  - 
which  motion  one  may  wish:to  consider  as  a special  type  of  lattice  im- 
perfection - can  be  treated  along  the  same  general  lines  as  the  magnetic 
dipole  interaction  in  these  cases.  Most  of  the  experimental  work  concerned 
with  this  type  of  thermal  motion  has  been  done  for  protons,  for  which  the 
quadrupole  interaction  is  of  course  absent.  Diffusion  in  solids  and  quadru- 
pole relaxation  has  already  been  discussed  earlier  in  this  review. 


F.  Excitons 

No  information  of  either  theoretical  or  experimental  nature  of  the 
quadrupole  interaction  around  these  imperfections  is  known  to  the  author. 
It  is  doubtful  whether  they  can  be  produced  in  sufficient  concentration  to 
cause  observable  effects. 


This  survey  may  be  terminated  with  these  concluding  remarks.  The 
influence  of  imperfections  on  nuclear  magnetic  resonance  signals  is  fre- 
quently so  profound  that  inversely  the  magnetic  resonance  may  be  used  as  a 
tool  for  the  investigation  of  the  nature  of  the  imperfections.  In  general,  the 
quadrupole  effects  are  more  pronounced  than  the  magnetic  effects.  Problems 
which  can  be  attacked  successfully  with  nuclear  resonance  techniques  include: 
ael^liilualori,  ionic  diffusion,  electronic  structure  of  « and  alloys, 

density  of  bound  and  free  carriers  in  semiconductors,  magnetic  properties 
of  F-centers,  density  and  distribution 'of  dislocations,  density  and  distribu- 
tion of  vacancies  and  interstitials,  lattice  deformation  and  electron  distribu- 
tion around  chemical  impurities,  short-range  order,  precipitation  and 
radiation  damage.  The  application  of  nuclear  resonance  techniques  to  these 
fields  has  only  just  started.  This  review  should  therefore  be  considered  as  a 
progress  report.  Many  new  results  are  likely  to  be  obtained  in  the  near 
future . 
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Appendix 

It  has  been  shown  early  in  the  history  of  nuclear  magnetic  relaxation 
[30}  that  the  relaxation  of  a spin  system  with  21  + 1 equidistant  levels  haB  the 
following  characteristics. 

The  saturation  under  the  influence  of  a ates-dy- state  radio-frequency 

£ m t 4 + m m - w—  mm 

X o ^ivcu  uy 

X"/X^  = An/Ano  = [ 1 y2H2f  Tlg(v  )] 

where  An  is  the  difference  in  population  between  adjacent  levels  in  the 
presence  of  an  oscillating  field  with  rms  amplitude  Hrf  and  AnQ  the  dif- 
ference under  thermal  equilibrium  in  the  absence  of  the  if  field.  The  re- 
covery of  the  sydtem  to  the  state  of  thermal  equilibrium  after  the  rf  field 
has  been  turned  off  is  given  by  a simple  exponential  function  with  character- 
istic time  T^. 

Whenever  the  splitting  of  the  spin  levels  is  non-equidiSant  the  problem 
of  relaxation  becomes  rather  complicated.  The  steady-state  condition  can 
still  be  analysed  readily  in  tegrms  of  th.'?  population  number  of  each  leyelj  if 
the  establishment  of  spin-spin  equilibrium  in  each  level  is  more  rapid  than 
the  spin-lattice  relaxation.  Pound  [19]  has  discussed  this  situation  for 
systems  with  quadrupole  coupling  and  Pake  [39]  for  the  case  of  magnetic 
hyperfine  splittings.  The  transitions  between  tin  different  energy  levels 
may  have  different  saturation  parameters.  The  recovery  is  now  a rather 
complicated  process  and  is  not  related  to  the  steady -state  saturation  in  the 
simple  manner  as  it  was  for  equidistant  levels.  In  the  case  of  quadrupole 
interaction  in  metals  the  recovery  can  be  analysed  in  detail.  We  shall 
describe  the  situation  which  applies  to  cold -worked  copper  (I  = 
external  field  produces  transitions  only  between  the  levels  irtj  = 

The  relaxation  process  is  entirely  due  to  interaction  with  the  conduction 
electrons  which  produces  transitions  between  adjacent  nuclear  spin  levels 
only.  The  transition  probability  per  unit  time  between  the  levels  m and  m - 1 
is 

w = Zn  r ( v ) U + m)(I  - m + 1 ) = { 2T  i ) ^(1  +m)(I  - m + 1) 


4).  The 
1 and  - T 
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It  is  independent  of  the  qn»drupole  splitting,  as  the  conduction  electrons 
take  up  any  small  energy  balance  in  their  kinetic  energy;  p^(v)  is  constant 
up  to  very  high  frequencies.  is  the  relaxation  time  in  well-annealed 
copper  where  the  levels  are  equidistant.  The  differential  equationsde scribing 
the  relaxation  process  are 


dn, 


3/2 


dt 


= + (-n,  /,  + n 


3/2  T "3/2  ' nl/2  ' “1/2 


° \ 
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dnl/2  .. 
dt 
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■at + x1  ~nl/2 
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Df  6 

+ n 


-1/2"  n - 1/2 
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° \ 
- n ,,,) 


Y2Hrf(_ill/2+n-l/2)«(v) 


The  third  equation  could  of  course  be  written  in  a form  symmetrical  to  the 
second,  but  it  shows  clearly  the  obvious  relation  that  the  total  number  of 
spins  n ^ + “i/z  + n 1/2  + n-3/2  remains  constant.  It  is  easily  verified  that 
under  steady -state  conditions  the  same  saturation  relation  is  still  satisfied. 
Except  for  the  email  change  in  v ) due  to  quadrupole  splitting,  the  nuclear 
signal  ih  cold-worked  copper  will  saturate  at  the  same  power  level  as  in 
annealed  copper. 

The  recovery  time  sfter  saturation  r»f  the  rentral  transition  is  different. 
For  H ^ = 0 the  system  of  equations  reduces  to  two  independent  linear  differen- 
tial equations  for  n^  - n ant*  n3/2  " n 3/  2 8°luti°r- 

nl/2  " n_  1/2  = Aexp{-t/Tl)  + Bexpt-bt/Tj)  + n®^  - n°y2 

The  constants  A and  B are  determined  from  the  intial  conditions  at  t = 0, 

n 1/2  ” n - 1/2  =0  *ndn3/2_n-3/2=n3/2'n-3/2‘n3/2  ^ ^ - 1/2  °r  dT  (n3/2_n-3/2)  = 0 
The  result  is  then  the  equation  used  in  the  text, 
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An  = | - | exp  (-t/Tj)  - | exp  (-6t/Tj)  + lj  AnQ 

Pulse  techniques  which  measure  thebe  transient  effect?  should  therefore  give 
different  results  in  the  coid-work’d  material  from  steady-state  techniques, 
unless  the  pulses  are  of  such  short  duration  and  consequently  their  frequency 
spectra  so  wide  that  they  turn  all  spins  rather  than  change  the  magnetization 
of  the  central  states  m^  = -j  or  - -y  only. 

If  the  dominant  relaxation  mechanism  is  caused  by  dipolar  interaction 
between  pairs  of  nuclear  dipoles,  the  set  of  differential  equations  becomes 
non-linear  and  would  not  permit  a solution  in  closed  form.  It  can,  how- 
ever, be  linearized  by  assuming  that  the  population  of  any  of  the  21  + 1 levels 

never  deviates  much  from  N/(2I  + 1).  Introduce  the  deviation  An  and 

m 

neglect  products  of  these  small  quantities.  In  the  case  of  quadrupole  relax- 
ation the  system  remains  linear  in  the  population  nuiiibcrs  juL  ouu  a solution 
for  transients  can  be  obtained  in  principle.  The  relation  to  the  steady-state 
solution  will  be  different  in  each  individual  ctise. 
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